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ABSTRACT 

The  conventional  low-power  CW  doppler  projectile  tracking  radar  is  investigated 
through  simulation.  The  power  requirement  for  full  range  tracking  with  real-time  track 
data  processing  stopped  within  160  m  of  the  point  of  impact  is  obtained.  This  power 
can  be  lowered  bv  a  factor  of  more  than  4  if  a  set  of  three  switchable  filters  is  utilized. 
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I.     INTRODUCTION 

At  a  proving  ground,  it  is  always  desirable  to  record  the  trajectories  of  test-fired 
munitions  for  later  analysis.  This  is  especially  true  when  anomalies  occur.  Continuous 
wave  (CW)  radars  have  often  been  adopted  for  this  purpose  because  of  their  simplicity. 
Compared  to  a  pulsed  radar  of  similar  range  performance,  a  CW  radar  has  no  modulator 
in  its  transmitter,  hence  no  range  gating  is  necessary  in  its  receiver.  Furthermore,  a  CW 
radar  has  a  bandwidth  corresponding  to  the  doppler  shift,  which  rarely  exceeds  100 
kilo-Hertz  (kHz)  while  a  pulsed  radar  has  a  band  width  in  the  mega-Hertz  range.  Thus 
a  CW  radar  can  be  designed  to  carry  out  coherent  signal  integration  over  a  longer  time 
interval.  This  increases  the  integration  gain  and  lowers  the  required  transmitter  power 
for  the  CW  radar. 

High  gain  antennas  of  very  narrow  beamwidths  are  also  used  to  reduce  the  trans- 
mitter power  requirement.  It  is  important  to  aim  such  antennas  accurately  on  the  target. 
To  acquire  the  bearing  and  elevation  angles  of  a  target  directly,  several  receiving  feeds 
together  with  complicated  comparator  and  multiple  data  channels  are  needed.  Since  the 
trajectory  of  a  test  round  can  be  roughly  predicted,  an  alternative  to  the  actual  acquisi- 
tion of  target  tracking  angles  is  by  steering  the  antenna  according  to  the  predicted  tra- 
jectory. The  directly  measurable  data  will  be  the  doppler  shift  of  the  target  as  a  function 
of  flight  time.  Discrepancies  between  this  data  and  the  doppler  shift  based  on  the  pre- 
dicted trajectory  should  be  carefully  analyzed  and  explained.  When  an  anomaly  occurs 
which  causes  the  target  to  break  off  track,  this  recorded  data  will  be  valuable  for  diag- 
nosing the  problem. 

Many  CW  radar  tracking  systems  are  available.  For  example,  a  model  DR  5S2  X- 
band  Instrumentation  Radar  manufactured  by  TERMA  Elektronik  AS  of  Denmark  is 
being  installed  at  the  Ta-Fu  Proving  Ground  in  Taiwan.  These  systems  typically  track 
a  long-range  projectile  only  up  to  the  summit  of  the  trajectory.  Since  the  point  of  impact 
of  a  long-range  projectile  is  of  great  interest,  it  is  desirable  to  extend  the  range  of  the 
radar  to  cover  as  large  a  portion  of  the  trajectory  as  possible  while  keeping  the  system 
relatively  simple.  The  purposes  of  this  thesis  are  to  investigate  the  maximum  range  at- 
tainable using  a  simple  CW  radar  and  to  determine  the  required  minimum  transmitter 
power. 


When  tracking  a  projectile,  a  CW  radar  detects  only  the  doppler  frequency  shift  due 
to  the  projectile  motion.  This  frequency  shift  is  proportional  to  the  component  of  the 
projectile  velocity  pointing  radially  away  from  the  radar.  By  combining  this  information 
with  the  tracking  angles,  the  ground  range  of  the  projectile  can  be  deduced  through  in- 
tegrating this  velocity  component.  A  low-pass  filter  can  be  employed  to  perform  this 
signal  integration  in  real  time.  The  output  can  be  used  to  update  the  predicted  track  if 
a  data  filtering  mechanism  such  as  the  Kalman  filter  is  utilized  in  real  time.  Since  in- 
formation on  the  doppler  integrating  filter  is  not  provided  by  any  manufacturer,  one  way 
to  implement  this  low  pass  filter  is  suggested  in  Chapter  2.  The  design  described  therein 
is  consistent  with  the  operational  instruction  for  the  aforementioned  system.  It  is  needed 
for  tracking  simulation  required  in  Chapter  IV.  The  trajectory-updating  Kalman  filter 
is  a  desirable  improvement  to  currently  available  systems  on  the  market.  Its  design  is 
outside  the  scope  of  this  thesis. 

As  explained  above,  trajectory  estimation  is  necessary  for  tracking  with  a  simple  CW 
radar.  Since  the  radar  and  the  gun  are  moving  with  the  spinning  Earth,  the  familiar 
reference  frame  in  which  the  radar  and  the  gun  are  stationary  is  an  accelerating  frame 
which  is  difficult  to  work  with.  Instead,  by  ignoring  the  rotation  of  the  Earth  around 
the  sun,  the  inertial  frame  of  reference  having  its  origin  at  the  center  of  the  Earth  and 
its  z-a\is  along  the  spin  axis  of  the  Earth  simplifies  the  equations  of  motion  of  the 
projectile.  The  z-x  plane  can  be  chosen  to  contain  the  radar  at  t  =  0  when  the  projectile 
is  fired.  To  obtain  the  familiar  ground-based  quantities,  transformation  to  the  inertial 
frame  in  which  the  radar  (or  the  gun)  is  passing  instantaneously  through  the  z-x  plane 
can  be  performed.  In  this  frame,  the  distance  and  the  relative  velocity  of  the  projectile 
from  the  radar  (or  the  gun)  can  be  found  easily. 

The  projectile  will  be  spinning  around  its  axis  if  a  rifled  gun  barrel  is  used.  The  lift 
force  due  to  air  pressure  will  rotate  the  projectile  to  align  its  axis  along  its  air  velocity 
vector.  The  kinematics  of  the  projectile,  including  its  angular  momentum,  is  most  easily 
described  in  terms  of  the  motion  of,  and  the  rotation  around,  its  center  of  mass  [Ref. 
1].  The  orientation  of  the  projectile,  which  lies  along  one  of  its  principal  axis  of  inertia, 
can  be  described  in  terms  of  two  of  the  three  Euler  angles  [Ref.  2],  leaving  the  third  angle 
for  measuring  the  spin  of  the  projectile. 

It  can  be  shown  that  the  effect  of  gravitational  force  on  the  projectile  orientation  is 
negligible.  The  ratio  between  the  gravitational  force  on  the  center  of  mass  to  the  force 
affecting  the  angular  momentum  is  of  the  order  of  the  square  of  the  ratio  of  projectile 
size  to  the  Earth  radius  [Ref.  3].   Thus  the  change  of  the  orientation  of  the  projectile  axis 


is  due  largely  to  the  resistance  of  air.  For  a  streamlined  body,  the  lift  from  the  air  flow 
is  large  when  there  is  an  angle  of  attack  [Rcf.  4].  Furthermore,  the  rotation  of  the 
projectile  axis  is  slow.  It  is  reasonable  to  assume  that  the  orientation  of  the  axis  of  the 
projectile  is  always  along  its  air  velocity.  This  assumption  will  be  utilized  to  determine 
the  radar  cross-section  of  the  projectile. 

It  is  easily  seen  in  the  ground  coordinate  that  the  projectile  velocity  is  much  greater 
than  the  wind  velocity  so  that  the  latter  only  introduces  small  corrections  to  the 
projectile  axis.  Since  the  wind  acts  on  the  projectile  over  the  entire  trajectory,  a  better 
estimate  of  the  range  can  be  obtained  by  including  the  wind  velocity  profile  to  trajectory 
prediction.  A  less  elaborate  alternative  is  to  adopt  a  wind  profile  model  such  as  the 
Ekman  spiral  for  wind  speed  approximation  [Ref  5].  Since  the  change  in  the  predicted 
range  is  not  expected  to  be  large,  the  wind  velocity  is  treated  as  constant  throughout  this 
thesis. 


II.     PROPOSED  C\\  RADAR  DESIGN 

A.     RADAR  SYSTEM  SCHEMATICS 

The  receiver  design  of  a  radar  is  critical  to  its  range  performance.  Since  no  specific 
information  for  any  CW  radar  tracking  system  is  available,  the  following  receiver  con- 
figuration is  proposed.      Here  Figure  1  shows  the  schematics  of  the  radar. 
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Figure   1. 


Schematics  for  Proposed  CW  Radar  Design. 


The  received  signal  at  the  frequency  (ft+fd)  from  the  receiving  antenna  is  fed  to  a 
narrow-band,  low-noise  RF  pre-amplifier  with  a  bandwidth  of  1  MHz  centered  at  10.475 
GHz.  The  noise  figure  of  such  a  pre-amplifier  can  be  as  low  as  1.5  to  2  dB,  with  a 
minimum  gain  of  28  dB  [  Ref  6  ].  This  low-noise  pre-amplifier  can  limit  the  receiver 
noise  figure  (  Fn  )  to  about  2  dB.  The  1  MHz  bandwidth  is  chosen  according  to  the  es- 
timated maximum  doppler  frequency  {fdm„)  and  the  fact  that  the  target  is  always  flying 
away  from  the  radar.  The  output  signal  from  the  low-noise  pre-amplifier  is  reduced  to 
the  base  band  at  the  mixer,  and  then  passed  through  the  video  filter  of  bandwidth  B,. 
The  filter  output,  at  the  doppler  frequency/,,  can  be  recorded  either  with  an  analog  or 
a  digital  recorder.  The  bandwidth  of  this  recorder  should  be  wider  than  fdmex  to  allow 
truthful  recording  of  the  raw  data  for  later  analysis.  The  doppler  signal  is  then  fed  into 
the  square-law  detector. 

The  output  of  the  square-law  detector  oscillates  at  a  frequency  of  2fd.  Since  this 
frequency  is  proportional  to  the  velocity  of  the  projectile  receding  away  from  the  radar, 
it  can  be  integrated  to  obtain  the  projectile  range  information. 

The  integration  is  most  easily  accomplished  by  counting  the  number  of  cross-over 
points  of  a  pre-determined  level  by  the  square-law  detector  output.  As  to  be  explained 
later,  this  number,  denoted  by  MX,  is  the  displacement  of  the  projectile  away  from  the 
radar  measured  in  half  the  transmitter  wavelengths. 

The  square-law  detector,  a  counter,  and  a  clock  will  be  adequate  for  measuring  the 
range  and  the  velocity  of  the  projectile  receding  from  the  radar.  First,  a  measurement 
base.  MB.  is  set  in  terms  of  the  number  MX  of  transmitter  half  wavelengths.  Then  the 
time  Tn  for  the  output  of  the  square-low  detector  to  run  through  MX  pre-set  level 

\jfJ3 

crossings  is  measured.    Here  MB  is  the  radial  range  increase  of  the  projectile  and  ' 
is  the  radial  velocity  of  the  projectile. 

As  discussed  in  the  introduction,  the  trajectory  of  the  projectile  is  predicted  before- 
hand. The  elevation  and  bearing  angles  of  the  radar  are  updated  at  the  end  of  every  Tn 
interval.  This  information  is  sent  to  the  servo  control  system  to  keep  the  antenna  on  the 
projectile. 

B.     THEORY 

The  theory  of  operation  for  the  proposed  CW  radar  system  will  be  described.  The 
doppler  effect  is  used  in  a  CW  radar  for  moving  target  detection.  This  effect  was  found 
by  Christian  Doppler  in  1842.  It  states  that  the  frequency  of  a  harmonic  signal  will  shift 
when  reflected  by  a  moving  target.    The  change  in  frequency  is  proportional  to  the  ve- 


locity  of  the  target.  Let  fd  denote  the  doppler  shift,  V,  the  radial  velocity  of  the  target 
with  respect  to  the  radar,  c  the  velocity  of  light,  and/  the  transmitted  carrier  frequency. 
The  doppler  frequency  shift  is 
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The  doppler  is  zero  when  the  trajectory  of  the  projectile  is  perpendicular  to  the  radar 
line  of  sight.    It  is  positive  for  an  approaching  target  (see  Figure  2). 
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Figure  2.      Spectra  of  Received  Signals:  (a)  no  doppler  shift,  no  relative  motion  (b) 
approaching  target  (c)  receding  target. 


The  frequency  of  the  signal  out  of  the  square-law  detector  is  2fd.    Thus  the  number 
of  crossings.  MX,  of  this  signal  over  a  pre-set  level  within   a  time  interval  Tn  is  : 


If 

MS  =f„x4i 
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Tnx2 
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According  to  equation  (2.1),  the  radial  displacement  of  the  projectile  over  this  period 
(Figure  3  on  page  7)  is  : 
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Figure  3.      Measurement  Base  and  Velocity  vs.  Time. 
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For  example,  if  the  transmitting  frequency/  of  interest  is  10525  MHz,  we  can  set 
the  number  MN  to  100.   Then 


MB  = 


100  x  2.997  x  W 
2  x  10525  x  106 


(2.7) 


=   1.4234  m 


(2.8) 


This  is  a  small  change  in  target  location  for  updating  the  radar  pointing  direction. 
Since  MB  is  the  displacement  of  the  target  along  the  radar  pointing  direction,  this  dis- 


placement  can  be  projected  onto  the  ground.   The  ground  range  can  be  derived  by  add- 
ing up  the  individual  projected  ranges. 


III.  PROJECTILE  TRAJECTORY  PREDICTION 

To  analyze  the  range  performance  of  the  radar,  information  about  the  trajectory  of 
the  projectile  is  required.  There  are  many  factors  such  as  air  drag,  gravitational  force, 
projectile  shape  and  spin,  Earth  rotation,. ..etc..  which  can  affect  the  trajectory.  As  dis- 
cussed in  the  INTRODUCTION,  the  projectile  spin  mainly  provides  a  stabilizing 
mechanism  to  the  projectile  and  is  assumed  to  have  little  effect  on  its  trajectory.  The 
direction  of  the  projectile  axis  (or  heading)  is  assumed  to  be  along  the  direction  of  the 
air  velocity  of  the  projectile.  Thus  only  the  motion  of  the  center  of  mass  of  the  projectile 
needs  to  be  considered.  The  projectile  center  of  mass  and  the  center  of  air  pressure  are 
assumed  to  be  the  same.   The  trajectory  is  thus  reduced  to  that  of  a  particle  trajectory. 

Wind  velocity  can  increase,  decrease  or  deviate  the  trajectory  through  air  drag  on 
the  projectile.  In  this  thesis,  it  is  assumed  that  wind  velocity  is  constant  over  the  course 
of  the  projectile;  it  is  usually  called  the  ballistic  wind. 

Drag  force  depends  on  the  air  velocity  of  the  projectile.  It  points  against  the 
projectile  heading.  Its  magnitude  is  proportional  to  the  square  of  projectile  air  speed, 
the  local  air  density,  the  projectile  cross-section,  and  the  drag  coefficient.  Air  density 
depends  on  the  temperature  profile  of  the  atmosphere.  The  drag  coefficient  is  very 
complicated  and.  in  fact,  varies  with  the  projectile  air  speed.  Measured  data  of  the  drag 
coefficient  are  available.   The  drag  force  on  the  in-flight  projectile  is  given  below  : 

f  =    -hj-ACD(V-Hf  (3.1) 

where  W  is  wind  velocity.  I 'is  the  velocity  of  the  projectile.  CD  is  the  drag  coefficient, 
A  is  the  cross-section  area  of  the  projectile,  and  p  is  the  local  air  density.  Here  h  is  the 
heading  of  the  projectile,  which  is  a  unit  vector  in  the  direction  of  V  —  IV, 

The  gravitational  force  is  the  only  other  factor  affecting  the  trajectory  in  an  inertial 
frame  of  reference.  In  the  accelerating  reference  frame  fixed  on  the  surface  of  the  Earth, 
there  is  the  fictitious  Coriolis  force.  It  causes  an  apparent  trajectory  drift.  As  observed 
on  the  ground,  this  drift  can  reach  3.1  km  for  a  trajectory  of  27.4  km  ground  range. 


A.     THE  TRAJECTORY  IN  A  FIXED  INERTIAL  FRAME 

In  order  to  simplify  the  equation  of  motion,  the  Earth-centered  inertial  frame  (r, 
8,  4>)  is  used.  The  plane  containing  the  radar  and  the  spin  axis  of  the  Earth  when  the 
projectile  is  fired  at  t  =  0  is  the  z-x  plane.  The  radar  is  moving  with  the  Earth  at  an 
angular  velocity  of  Q.  .    Its  location  in  this  inertial  frame  is  (rR,  8R,  Q/)  . 

The  components  of  the  projectile  velocity  in  this  reference  frame  are  : 

r=Vr  (3.2) 

r8=Ve  (3.3) 

r(  sin  0)0=^  (3-4) 

The  Lagrangian  formulation  of  the  principle  of  least  action  is  convenient  for  deriv- 
ing the  equations  of  motion  of  the  projectile  [Ref.  7].  The  Lagrangian  of  the  projectile 
in  the  inertial  spherical  coordinate  system  [Ref.  1  ,  8]  (  r,  8,  </>  ),  without  considering  the 
drag  force,  can  be  expressed  as  : 

L=  T{r,6,<j>,r,e,4>)-  U{r,6,4>)  (3.5) 

L  =  -y-  (P  +  r262  +  r2<p2  s\r?8)  -  mga{  1  -  -f )  (3.6) 

where  g  is  the  gravitational  acceleration  at  ground  level.     For  the  r  component,  the 
equation  of  motion  is  : 

d    cL  cL  (3?) 


di     cr  cr 


1 

r=    --*y  g  +  re2+rj>2  sin2 6  (3.8) 

r 


For  the  8  component  : 

d    cL  cL 

di    £@  66 


(3.9) 


-^(rtf)  =  r2i>2  sin  0  cos  8  (3.10) 

di 


10 


and 

,  for  the 

0  component 

d 

cL 

cL 

dt 

4> 

2 

C(p 

(3.11) 


-^-(-^c/>sm20)  =  O  (3.12) 

The  nonlinear  drag  force  can  be  included.    The  complete  equations  of  motion  become  : 

I  '2  4-    V2 
dl 


K=  -(-frs  +  — j — +fr  (3-13) 


-^(■7Tve)  =    ™tdvl  +  rfe  (3.14) 

-^-(rsmdVJ  =  rsin6|/i  (3.15) 

r  =  o  +  /j  (3.16) 

where  h  is  the  height  from  the  Earth  surface  and  a  is  the  radius  of  the  Earth.    The  ratio 
of  projectile  height  to  the  Earth  radius  can  be  expressed  as  : 

T   =  1+4  0-17) 
The  projectile  velocity  is  : 

V  =    Vrr+  Ved+  VJb  (3.18) 

=  rr  +  rd6  +  r  sin  6(j>4)  (3.19) 
and  the  wind  velocity  is  : 

W  =   lVe6+  wj<j>+  WTr  (3.20) 

Note  that  fr,  fe,  fc  are  the  components  of  the  drag  force  divided  by  the  mass  of  the 
projectile  : 


/  =    -h-f-ACD\V-  IV  \2  (3.21) 

2m         u  v 
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-  ff  +fe6  +fo<t>  (3-22) 

The  equations  of  motion  can  be  discretized.  Assume  that  the  gravitational  and  the 
drag  force  acceleration  are  approximately  constant  over  the  period  T„  when  the  projectile 
moves  from  r„  to  »•„_,.  the  equations  of  motion  are  approximately  : 

=    -g£-.  +  rnen2  +  rnsin2dri<!>l+frn  (3.23) 


Tn 


2  2  ' 

rn+}6n+]  -  rn$  2  . 2 

—    =  r„  sin  0n  cos  $„<<>„  +  r,fen  (3. 24] 


T 

1  r. 

r2n+]  sm2dn+l4>n+]  -  r2n  sm2en4>„ 


=  rnsm0nLn  (3.25) 


1  n 

From  the  above  equations,  if  r„,  QR,  (/>„  and  r„,  8„.  4>n  are  known,  then  r„.u  0nJ.j,  (/>„_,. 
rn+u  6n_u  and  0„_,  can  be  derived  as  follows  : 

With  measured  Tr,  equation  (1)  provides  >■„_,  »  then 

rn+]  =  /,  + 2 r"  (3,26) 

Equations  (3.26)  and  (3.24)  give  8„^u  then 

en+]  =  en  +  — j^-Tn  (3.27) 

Equations  (3.27)  and  (3.23)  give  <£„_,.  then 

^■H    =    <^n+  2  T"  (3-28) 

Note  that  the  muzzle  location  is  ~rc  and  the  muzzle  velocity  of  the  projectile,  which 
is  assumed  to  be  known,  is  Vc.  Through  the  above  procedures,  the  full  trajectory  of  the 
projectile  can  be  constructed  by  computer  simulation. 
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B.     GROUND-OBSERVED  TRAJECTORY 

In  order  to  obtain  the  trajectory  as  observed  from  the  ground,  the  position  of  the 
projectile  should  be  converted  into  its  relative  location  from  the  radar  antenna  and  from 
the  gun.  The  necessary  coordinate  transformation  matrices  are  given  below  (Figure  4 
on  page  13). 


Figure  4.      Inertial  Coordinate  Systems. 


1.     Trajectory  Relative  to  Radar 

To  set  the  elevation  and  bearing  angles  for  the  radar  and  to  estimate  the  doppler 
shift,  the  location  and  velocity  of  the  projectile  relative  to  the  radar  as  viewed  from  the 
ground  are  needed.  On  the  ground,  the  projectile  motion  is  most  easily  described  in  an 
inertial  frame  in  which,  at  time,  t.  the  radar  is  passing  through  the  z  —  x'  plane  instan- 
taneously. The  z'  axis  is  chosen  so  that  it  coincides  with  the  z  axis.  The  transformation 
from  the  fixed  inertial  spherical  coordinate  system  to  this  (x',y'  ,z')  system  can  be  ac- 
complished through  the  following  rotation  matrices  R„  Rrec  ,  and  Rcy,  : 
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RJLQt)  = 


cos  ilt      sin  Qt     0 

-  sin  Qi     cos  Q/     0 

0  0         1 


(3-29) 


RrecW   = 


cos  4>p     0     -  sin  ^ 

sin  <t>p     0      cos  0f 

0  1  0 


(3.30) 


Rcytfr)  = 


sin  ^       cos  flr      0 

cos  8p     —  sin  8p     0 


0 


0  1 


(3.31) 


a.      Velocity  Transformation 

The  projectile  velocity  in  the  (x',y',  z')  coordinate  system  can  be  obtained 
by  first  changing  the  results  in  section  II1A  from  the  spherical  coordinate  system  into 
the  rectancular  coordinate  svstem  (x,  v.  z)  : 


px 


py 
VP2 


=  Rrec(<t>p)Rcy^r] 


I' 


Po 


! 


pc 


(3.32) 


Since  the  (x',y',  z')  coordinate  system  is  obtained  by  rotating  the  (x,  y,  z) 
coordinate  svstem  bv  an  ansle  Qi  around  the  z-axis. 


V. 


px 


V, 


py 


\\ 


p* 


=  RJQi 


V. 


f  '■ 


V, 


py 


p: 


(3-33) 


Rrec(d>p  -  Qt)Rcy£6p) 


pr 


p6 


r 


PC 


(3.34) 


The  velocity  of  the  radar,  \'R,  in  the  (x',  y'  ,z')  coordinate  is 
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=  RrK{O)RU0R 


0 
0 


(3.35) 


sin  6R      cos  0R      0 

0  0  1 

cos  6R     —  sin  0R     0 


0 
0 


(3.36) 


0 
0 


(3.37) 


Hence 


VD  =  r 


*;'•> 


=  Q.{rRsmBR)y- 


b.     Projectile  Location  Relative  to  the  Radar 

The  relationship  between  (x',y',  z')  and  (r,  6,  <f>)  is  given  by 


(3.38) 
(3.39) 


X 

r  sin  6 

COS  <f) 

V 

= 

r  sin  6  sin  0 

-2  - 

r  cos  0 

~x'~ 

X 

t 

V 

t 

2 

=  RJfr) 

y 

r  sin  6  cos(</>  —  Qt) 

= 

r  sin  6  sin(0  —  Q/) 

r  CC 

sd 

(3.40) 


(3.41) 


(3.42) 
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c.     Radar  Coordinate 

Substitute  (rR,  9R.  Qt)  into  equation  (3.42), 


*R 

rR  sin  6R 

>'r 

= 

0 

ZR 

rR  cos  0R 

(3.43) 


2.     Trajectory  Relative  to  Gun 

Relative  to  the  gun  (muzzle),  a  similar  inertial  frame  (x  "  .y  "  ,  z  "  )  in  which  the 
gun  is  instantaneously  passing  through  the  z  "  -  x  "  plane  at  time  t  can  be  used.  The 
z  "   axis  coincides  with  the  z-axis.    The  coordinate  transformation  is  : 


;•  sin  6  cos(0  -  4>G  —  Q.I) 

r  sin  6  sin((/>  —  4>G  —  Qi) 

r  cos  6 


(3.44) 


The  firing  plane  of  the  gun  is  the  plane  containing  its  barrel  and  passing  through 
the  center  of  the  Earth.  Tor  artillery,  it  is  customary  to  use  the  coordinate  system 
(xf,yf,zf)  in  which  the  gun  muzzle  is  the  origin  and  the  zf—  xf  plane  is  the  firing  plane. 
Here  z,  is  the  vertical  direction  and  x,  is  along  the  horizontal  from  the  gun.  The  relation 
between  (xf.  v>,  zf)  and  (x  "  ,y  "  ,z  ")  depends  on  the  gun  bearing  angle  8A!G.  which  is  the 
angle  of  the  x,  axis  measured  counterclockwise  from  the  East.  If  the  the  gun  is  oriented 
along  8AZC  and  is  located  at  rc  ,  60.  <f>a,  the  vector  r  —  rc  in  the  coordinate  (xf,y„  z,)  is  : 


0      cos  dAzG       sin  BA:G 

0  -  sin  0AzC     cos  6AzC 

1  0  0 


cos  6, 


0     sin  0 


1 


0 


-  sin  6G     0     cos  0 


- 

.-            -J 

0 

X  " 

y  " 

c  _ 

(3.45) 


Note  that  /.  —  rG  is  the  relative  location  of  the  projectile  from  the  gun.    Hence 
its  components  can  replace  (x  "  ,y     ,  z  "  )  direct!}'  in  equation  (3.45). 
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IV.     RANGE  ANALYSIS  AND  POWER  REQUIREMENT 

A.     RANGE  EQUATION 

The  power  requirement  of  the  CW  radar  can  be  obtained  from  the  radar  equation 
[Ref.  9]  : 

(4n)3kTsFnL(-$r)min(BvR4) 

p.  =  ~rh C4.I) 

G  oA  t 

In  equation  (4.1),  the  radar  cross-section  {a),  the  maximum  range  (R),  and  the 
propagation  factor  (F)  will  be  discussed  in  this  chapter  in  detail.  The  Boltzman  constant 
k  is  287.05  n-m  kg-°A'.  The  absolute  temperature  Ts  is  290°A'  .  The  wavelength  A  is 
0.0285m.  Other  parameters  are  based  on  the  configuration  of  the  CW  radar  as  proposed 
in  Chapter  2.   They  are  defined  as  follows. 

1.  The  Receiver  Noise  Figure 

The  receiver  noise  figure  (Fn)  is  the  ratio  of  the  signal-to-noise  power  ratio  at  the 
receiver  input  to  the  signal-to-noise  power  ratio  at  the  receiver  output,  assuming  that  the 
receiver  temperature  is  290eA'  .  The  Fr  value  for  the  receiver  considered  is  assumed  to 
be  2  dB. 

2.  System   Loss   Factor 

The  loss  factor  consists  of  the  fractional  loss  of  transmitter  power  delivered  to 
the  transmitting  antenna  (A,),  the  fractional  loss  of  the  receiving  antenna  output  power 
delivered  to  the  receiver  input  (A,),  and  miscellaneous  loss  (Lm).  They  are  discussed  be- 
low. 

a.      Transmission  Loss 

The  transmission  loss  is  the  loss  in  the  transmitting  or  the  receiving  portion 
of  the  transmission  line.  Since  the  transmitter  is  light-weight  and  is  to  be  mounted  very- 
close  to  the  antenna,  L,  is  estimated  to  be  about  1.2  dB  due  to  the  use  of  narrow  band- 
width rotary  joints.  The  receiver  pre-amplifier  can  be  mounted  directly  to  the  receiving 
antenna  output,  followed  with  the  mixer  immediately.  This  arrangement  will  reduce  A, 
to  about  the  same  level  as  A,. 
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b.     Miscellaneous  Loss 

The  miscellaneous  loss  (Lm)  for  this  C\V  radar  system  is  mainly  due  to  the 
quantization  loss  at  the  comparator.  Since  only  two-level  (crossing)  comparison  is 
made,  it  is  not  expected  to  be  significant  and  will  be  ignored. 

Base  on  the  above  argument,  the  total  loss  is  estimated  at  2.4  dB. 
3.     Signal-to-Noise  Ratio 

For  a  tracking  radar,  it  is  important  that  the  receiver  output  signal  is  sufficiently 
high  above  the  noise  level  so  that  information  carried  by  the  signal  can  be  extracted  re- 
liably. With  the  proposed  receiver  configuration,  the  signal-to-noise  ratio  (S'N)  of  the 
output  from  the  video  filter  is  the  figure  of  significance.  First  of  all,  this  signal  is  re- 
corded for  later  analysis.  Its  minimum  S  N  has  to  be  assured.  In  real-time,  after  passing 
through  the  square-law  detector,  this  signal  generates  the  timing  sequence  for  updating 
the  radar  pointing  direction  and  for  evaluating  the  projectile  radial  velocity.  It  drives  the 
gain  control  to  keep  the  reference  signal  at  a  preset  level.  It  is  also  the  input  into  the 
comparator  following  which  the  number  of  signal  cross-overs  is  tallied  and  the  elapsed 
time  is  clocked  as  explained  in  section  1 1  A. 

The  square-law  detector  and  the  comparator  perform  threshold  detection.  The 
reference  signal  acts  as  the  threshold.  The  comparator  outputs  a  +  1  signal  when  the 
detector  output  is  above  the  reference  level:  it  outputs  a  -1  signal  when  the  detector 
output  falls  below  the  reference.  The  reference  signal  level,  which  is  set  at  a  fraction  of 
the  detector  output  signal,  determines  the  accuracy  of  the  real-time  output  of  the  system. 
In  what  follows,  the  minimum  S  N  ratio  and  the  threshold  setting  in  terms  of  the 
threshold-to-signal  power  ratio  are  discussed. 

The  signal  input  to  the  comparator  is  oscillating  at  twice  the  doppler  frequency. 
Within  each  cycle  of  the  signal,  the  portion  rising  above  the  threshold  level  to  the  peak 
and  back  to  the  threshold  level  can  be  considered  as  a  pulse.  This  pulse  is  contaminated 
with  noise,  which  is  assumed  to  be  white.  The  probability  of  each  contaminated  pulse 
falling  under  a  certain  threshold  level  is  well-known.  This  probability  of  missed  de- 
tection is  the  fraction  of  signal  pulses  failed  to  be  detected  by  the  comparator.  Between 
two  pulses,  the  number  of  false  alarms  is  the  number  of  extra  cycles  counted.  These  two 
effects  tend  to  cancel  each  other  so  the  optimum  threshold  setting  is  achieved  by  equal- 
izing the  fraction  of  missed  pulse  and  the  fraction  of  noise  pulses.  The  program  to 
compute  these  fractions  are  listed  in  Appendix  D  [Ref.  10]. 

Figure  5  on  page  19,  Figure  6  on  page  20,  Figure  7  on  page  20.  Figure  8  on 
page  21,  and  Figure  9  on  page  21  show  the  fractional  errors  due  to  missed  pulses  and 
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noise  pulses  respectively.  If  the  signal-to-noise  ratio  is  10  dB,  the  comparator  output 
error  is  5.55%  when  the  threshold-to-signal  power  ratio  is  set  at  0.455.  Stronger  signal 
allows  the  threshold  to  be  set  at  a  lower  fractional  value  with  an  improved  error  ratio. 
For  example,  for  a  signal  12  dB  above  the  noise,  the  threshold  setting  of  0.38.  or  -4.2 
dB  below  the  signal,  will  limit  the  error  to  1.15%.  In  this  thesis,  the  minimum  signal- 
to-noise  ratio  is  chosen  to  be  10  dB. 
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Figure  5.      Missed  Detection  and  Noise  Error  vs.  Threshold-to-Signal  Poner  Ratio 
S/N-  10  dB. 
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Figure  6.      Missed  Detection  and  Noise  Error  vs.  Threshold  to  Signal  Power  Ratio 


SN=   10.5  dB. 
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Figure  7.      Missed  Detection  and  Noise  Error  vs.  Threshold  to  Signal  Power  Ratio 
SN=   11  dB. 
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Figure  8.      Missed  Detection  and  Noise  Error  vs.  Threshold-to-Signal  Power  Ratio 
SN=  11.5  dB. 
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4.  Transmitting  and  Receiving  Antenna  Gain 

Separate  hut  identical  transmitting  and  receiving  antennas  are  used  to  reduce 
direct  coupling  of  the  transmitter  power  into  the  receiver.  Since  this  separation  is  small, 
it  is  neglected  in  the  analysis.  In  x-band.  antennas  with  one-degree  beamwidth  are  easily 
achievable.    Based  on  this  beamwidth,  the  antenna  gain  (G)  [Ref.  9]  is  : 

c  =  -22222.  (4.2) 

=  AA.MiB  (4.3) 

and  the  diameter  D  of  such  an  antenna  is 

D  - 72  *  is  (4-4j 

=  72  x  M^l  (4  5) 

=  2m  (4.6) 

5.  Video  Bandwidth  and  Filters 

The  required  minimum  video  bandwidth  (Z?J  is  a  function  of  the  doppler  fre- 
quency. As  previous  mentioned,  the  bandwidth  should  be  greater  \hanfd„,ox  •  Although 
the  doppler  frequency/,  of  the  projectile  increases  rapidly  right  after  firing,  it  decreases 
monotonically  afterwards  by  a  factor  of  more  than  4  over  its  trajectory.  An  example  is 
shown  in  Figure  10  on  page  23. 
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Figure   10.      Doppler  Frequency  vs.  Slant  Range. 

Referring  to  the  range  equation  (-4.1),  in  order  to  use  the  lowest  transmitted 
power  possible  for  full  range  tracking,  the  radar  range  and  video  bandwidth  can  be 
traded  off.  Note  that,  assuming  a  constant  radar  cross-section,  the  transmitter  power 
depends  on  the  product  R'B^.  A  filter  of  narrower  bandwidth  may  be  employed  when 
the  projectile  reaches  a  greater  distance.  Adoption  of  a  set  of  three  switchable  doppler 
filters  is  recommended.  The  following  is  an  example  for  its  implementation  based  on  the 
doppler  frequency  of  the  projectile  given  in  Figure  10. 

The  maximum  doppler  shift  fdmox  from  the  projectile  shown  in  Figure  10  is  a 
little  less  than  54  kHz.  The  maximum  range  R  is  about  27  km.  Three  filters  are  to  be 
used  over  three  different  ranges.  Each  of  the  filters  will  have  a  bandwidth  of  1.2  times 
the  maximum  expected  doppler  frequency  over  the  range  in  which  it  is  to  be  used  to  al- 
low for  a  20%  estimation  error.  The  transmitter  power  level  will  be  minimized  if  the 
fourth  power  of  the  maximum  range  times  the  maximum  doppler  shift  in  the  corre- 
sponding range  is  a  constant. 
a.     First  Filter 

When  the  projectile  is  within  16.06  km  slant  range,  the  video  bandwidth  is 
chosen  to  be  1.2  x.fdmBx.   Thus  : 
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Bv  =   1.2  x^  (4.7) 

=  63.9  kHz  (4.8) 

b.  Second  Filter 

When  the  projectile  moves  beyond  16.06  km.  but  within  23.73  km  slant 
range,  the  doppler  frequency  has  a  maximum  value  of  about  22.03  kHz.  The  video 
bandwidth  is  : 

Bv  =  1.2  x  22.03  kHz  (4.9) 

=  26.44  kHz  (4.10) 

c .  Third  Filter 

When  the  projectile  reaches  23.73  km  slant  range,  the  maximum  doppler 
frequency  is  about  12.7  kHz.    The  video  bandwidth  is  : 

Bv  =  1.2  x  12.7  kHz  (4.11) 

=   15.2  kHz  (4.12) 

Since  RAfdmax  is  now  less  than  1  4  of  the  original  value  when  only  one  video  filter  is  used, 
the  minimum  required  transmitter  power  will  be  reduced  by  more  than  a  factor  of  4. 
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B.     TRAJECTORY  SIMULATION  AND  RESULTS 

Based  on  the  procedures  for  trajectory  estimation  presented  in  Chapter  III.  a  simu- 
lation program  was  written  and  is  included  as  Appendix  A.   It  was  written  in  PC-Matlah 
3S6.   The  program  is  described  as  follows. 
1.     Input  and  Parameters. 

a.  Earth 

The  operation  is  assumed  to  take  place  in  Taiwan.  The  radar  is  located  at 
the  longitude  of  121°36'27  "  East  and  the  latitude  of  24°36'3S  "  North.  The  Earth  ra- 
dius is  approximately  6.3677225  x  10'  m  (average  value  of  the  artic  radius 
6.356912  x  106  m  and  the  equator  radius  6.37S533  x  106  m).  The  gravitational  acceler- 
ation at  sea  level  is  9. SI  m/ sec-\  The  angular  velocity  of  Earth  spin  is 
7.2921  lS44x  10-'  rad  sec  [Ref.  11]. 

b.  Radar 

The  transmitting  frequency  is  10525  MHz.  The  measurement  base  MB  can 
be  varied  to  change  the  update  rate  of  the  antenna  pointing  direction.  For  both  the  re- 
ceiving and  the  transmitting  antenna,  the  beamwidth  is  lc  and  the  gain  is  44.3  dB.  The 
antenna  height  and  the  distance  between  the  gun  and  the  radar  can  affect  the  initial 
readings  of  the  doppler  signal. 

c.  Gun 

The  firing  elevation  angle  is  assumed  to  be  45°  and  the  azimuth  is  counter- 
clockwise from  the  East  and  is  also  assumed  to  be  45°.  Both  the  elevation  angle  and  the 
azimuth  angle  can  be  varied.  The  distance  between  the  gun  and  the  radar  can  be  ad- 
justed. (Variations  of  up  to  2000  m  are  possible  at  the  test  site.)  The  barrel  length  of 
the  gun  under  consideration  is  7  m. 

d.  Projectile 

The  target  is  a  155  mm  calibre  artillery  projectile.  Its  weight  is  45.4  kg,  with 
a  total  length  of  938  mm.   The  maximum  muzzle  velocity  is  850  msec. 

e.  Environment 

(1)  Temperature  and  air  density  profiles.  The  temperature  Tc  at  sea  sur- 
face is  assumed  to  be  28S°A'  and  the  temperature  at  the  projectile  height  h  is  based  on 
the  following  profile  [Ref.  12]  : 

Th  =   T0-ajh  (4.13) 
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where  h  is  the  height  from  sea  level,  aT  is  the  temperature  lapse  rate  (-0.0065  °K  m). 
The  temperature  (7j  is  used  for  calculating  the  air  density  (pj  at  the  projectile: 

p*  =  po(-t>c~(*,+,:  <^4> 

1  o 

In  equation  (4.14)  p0  is  the  air  density  at  sea  level  (1.225  kg  m3),  and  R  is  the  gas  con- 
stant (287.05  n-m  kg-°A~). 

(2 j  Wind  profile.  The  wind  is  assumed  to  be  constant  during  the  flight 
time  over  the  projectile  trajectory.  The  value  used  is  the  so-called  ballistic  wind.  It  can 
be  expressed  as  : 

W  =    llrr  +  lVe6  +  WJ>  (4.15) 

In  this  thesis,  the  wind  is  6  m  sec,  North-West.  It  should  be  noted  that  the  extent  to 
which  the  estimated  trajectory  can  be  modified  by  the  actual  wind  profile  is  subject  to 
further  study. 

(3)  Mach  number.  The  Mach  number  is  the  ratio  o[  the  projectile  air 
speed  to  local  sound  speed.    The  sound  speed  at  the  height  of  the  projectile  is  : 


ah  =  JKRTh  (4.16) 

where  K  is  1.4  for  idea  gas.    The  Mach  number  (M)  is  given  by  : 

I  V  -  IV  | 
M  = ; ~  (4.18) 

where  |  V—  W\  is  the  projectile  speed  with  respect  to  air. 
f.     Drag  Coefficient 

The  drag  coefficient  (Q,)  is  required  for  calculating  the  drag  force.  It  de- 
pends on  the  shape  of  the  projectile.  The  CD  curve  is  measured  experimentally.  It  is 
given  as  a  function  of  Mach  number  M.  For  a  155mm  projectile,  the  CD  curve  is  shown 
in  Figure  1 1  on  page  27. 

2.     Results  of  Simulation 

The  results  of  trajectory  simulation    are  described  below.    The  variable  names 
used  in  the  program  are  given  in  the  parentheses  when  they  are  first  mentioned. 
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a.  Radar  Observed  Trajectory 

The  slant  range  (SRR)  is  the  distance  from  the  radar  antenna  to  the  in-flight 
projectile.  The  curve  showing  the  projectile  height  (Height)  versus  its  slant  range  from 
the  radar  is  shown  in  Figure  12  on  page  2S.  The  total  flight  time  is  87.4  seconds,  and 
the  maximum  range  is  27.4  km. 

b.  Gun  Based  Trajectory 

The  ground  range  is  measured  from  the  gun  to  the  projectile.  The  projectile 
height  vs.  in-flight  time  (t)  and  the  projectile  height  vs.  ground  range  are  shown  in 
Figure  13  on  page  28  and  Figure  14  on  page  29.  The  maximum  ground  range  (Grange) 
is  27.4  km. 
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Drag  Coefficient  vs.  Mach  Number  for  155  mm  Projectile. 
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Figure   12.      Projectile  Trajectory:  Height  vs.  Slant  Range. 
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Figure   13.      Projectile  Trajectory:  Height  vs.  Time. 
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Figure  14.      Projectile  Trajectory:  Height  vs.  Ground  Range. 

The  cross  range  (CRange)  measures  the  drift  of  the  projectile  from  the  firing 
plane.  The  firing  plane  is  the  vertical  plane  containing  the  gun  barrel.  In  the  case  con- 
sidered, the  maximum  cross  range  is  about  6.3  km  South-East  of  the  firing  plane.  The 
point  of  impact  is  about  3.14  km  away  from  the  firing  plane.  The  curve  of  cross  range 
vs.  ground  range  is  shown  in  Figure  15  on  page  30. 
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Figure   15.      Projectile  Trajectory:  Cross  Range  vs.  Ground  Range. 

c.     Radar  Tracking  Parameters 

The  azimuth  (traz)  and  elevation  (trel)  antenna  pointing  angles  are  required 
for  tracking.  The  curves  of  these  angles  vs.  in-flight  time  when  the  gun-radar  distance 
is  500  m  are  shown  in  Figure  16  on  page  31  and  Figure  17  on  page  31. 

The  azimuth  and  elevation  angular  speeds  of  the  antenna  are  shown  in 
Figure   IS  on  page  32,  and  Figure   19  on  page  32. 
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Figure   17.      Radar  Elevation  Angle  vs.  Time. 
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The  antenna  servo  system  has  to  point  the  antenna  to  within  half  a  beamwidth  in  the 
direction  of  the  projectile.  Thus  the  half  beamwidth  divided  by  the  angular  speed  of  the 
antenna  is  the  required  response  time  of  the  servo  system.  The  servo  response  time  vs. 
the  projectile  in-flight  time  is  shown  in  Figure  20  on  page  33.  From  the  above  results, 
it  can  be  seen  that  over  the  first  few  seconds,  the  elevation  angular  speed  is  very  high. 
The  distance  between  the  gun  and  the  radar  can  be  increased  to  reduce  this  speed. 
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Figure  20.      Servo  System  Response  Time  vs.  Time. 

d.     Doppler  Frequency  and  Velocity 

The  doppler  frequency  (fd)  is  a  directly  measurable  data  by  a  CW  radar. 
The  doppler  cycle  can  be  detected  with  the  square  law  detector.  In  the  proposed  CW 
radar  configuration  of  this  thesis,  the  doppler  cycle  is  derived  from  the  time  interval  Tn 
for  each  measurement  base  traveled.  This  timing  sequence  has  to  be  simulated  in  the 
tracking  program.  In  the  simulation  program,  the  doppler  velocity  (Vdcpp)  is  calculated 
from  the  projectile  velocity  relative  to  the  radar  and  the  radar  pointing  direction.  Then 
Tn  is  estimated  from  the  time  required  to  travel  the  distance  MB  with  this  doppler  ve- 
locity. The  curves  of  doppler  frequency  vs.  slant  range,  doppler  frequency  vs.  time,  and 
doppler  velocity  vs.  time  are  shown  in  Figure  10  on  page  23,  Figure  22  on  page  34,  and 
Figure  21  on  page  34. 
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Based  on  the  above  results  from  trajectory  simulation,  the  projectile  radar 
cross  section  and  the  transmitter  power  can  be  estimated.   They  will  be  discussed  later. 


Figure  21.      Doppler  Frequency  vs.  Time. 
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Figure  22.      Doppler  Velocity  vs.  Time. 
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C.     PROJECTILE  RADAR  CROSS-SECTION 

The  projectile  radar  cross-section  (RCS)  is  denoted  by  o.  It  is  the  area  intercepting 
the  amount  of  plane  wave  power  which,  when  scattered  equally  in  all  directions, 
produces  an  echo  at  the  radar.    In  other  words,  the  RCS  is  given  by  : 

£ 

a  =  A7tr2{  -p-  f  (4.  IS) 

where  Es  is  the  scattered  electric  field  strength  from  the  target,  r  is  the  distance  from  the 
target  in  its  far  field  where  the  quantity  Es  is  measured,  and  E,  is  the  incident  plane  wave 
field  strength  at  the  target. 

Determination  of  the  RCS  of  a  complex  target  is  difficult.  Usually  the  radar  cross- 
section  of  a  target  depends  on  the  following  factors:  a)  target  size,  b)  target  shape,  c) 
radar  frequency,  d)  target  aspect  angle,  and  e)  radar  polarization.  The  RCS  is  greatly 
simplified  if  the  target  size  is  large  compared  to  the  radar  wavelength.  Consider  the  ra- 
dar cross-section  of  a  simple  sphere.  The  RCS  as  a  function  of  the  circumference  of  the 
sphere  measured  in  wavelengths  {2naj}.)  is  shown  in  Figure  23  on  page  36.  The  region 
where  the  target  size  is  large  compared  to  the  radar  wavelength,  i.e..  where  2naj).  is 
greater  than  10,  is  called  the  optical  region  [Ref.  9].  The  RCS  in  this  region  is  roughly 
the  physical  cross-section  at  area  of  the  target. 

The  155  mm  projectile  is  in  the  optical  region  at  the  transmitter  frequency 
(10525  MHz)  of  interest.  Its  RCS  value  will  be  approximated  with  its  geometrical 
cross-sectional  area.  Because  of  the  rotational  symmetry  of  the  projectile,  the  radar 
cross-section  of  the  projectile  is  a  function  of  the  angle  6asp  between  the  radar  pointing 
direction  {D)  and  the  projectile  heading  (/;)• 

The  radar  pointing  direction  D  is  given  by  : 


D=     _J       S  (4.19) 

where  rP  is  the  projectile  position  vector,  rR  is  the  radar  position  vector  as  defined  in 
Chapter  III. 
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The  heading  of  the  projectile,  h,  is  in  the  direction  of  its  air  velocity.    It  is  also  de- 
fined in  Chapterlll: 


(V-  U\) 
h  =— ~- 

W-  »'  I 


(4.20) 


In  terms  of  h  and  D,  the  aspect  angle  6o;p  of  the  projectile  as  seen  from  the  radar  is 


asp 


=  cos  \D.h) 


(4.21) 


Note  that  0  <  6gsp  <  n.  For  simplicity,  the  projectile  is  considered  as  a  circular  cone  on 
top  of  a  circular  cylinder,  as  shown  in  Figure  24  on  page  37.  The  height  of  the  cylinder 
is  414.4  mm.  The  height  of  the  cone  is  523.6  mm.  The  diameter  of  the  cylinder  and  the 
base  of  the  cone  is  155  mm. 
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Figure  24.      Approximate  Shape  of  155mm  Projectile. 


The  vertex  angle  of  the  cone,  2a,  is  given  by: 


o.  =  tan 


-l 


155 


2  x  523.6 


=  0.1469  rod. 

=  8.42° 


(4.22) 

(4.23) 

(4.24) 


The  physical  cross-section  of  the  cone  is  the  projected  area  of  the  projectile  along  D.  the 
radar  pointing  direction,  onto  a  plane  perpendicular  to  D.   This  cross-sectional  area  will 
be  discussed  below  according  to  whether  |  D  •  h  \  <  cos  y.  or  |  D  •  /?  |  >  cos  a. 
1.     Case  1 

As    \  D  •  h\>  cos  a,   the   projection    of  the   projectile   along   D   is    shown   in 
Figure  25  on  page  38. 
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Figure  25.      Radar  Cross  Section  Case  1. 


The  total  RCS  value  is 


<7,  =  na~  \  D  .  /?  |  +  lab  |  D  x  /j 


(4.25) 


where  <7  in  the  radius  of  the  cylinder  and  b  is  the  height  of  the  cylinder. 
2.     Case  2 

As  \D»h\  <  cos  a,  the  total  RCS  is  given  by  equation  (4.20)  (Figure  26  on 
page  39). 
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Figure  26.      Radar  Cross  Section  Case  2. 


o,  =  CT]  4-  a2  +  o3  +  oA 


(4.26; 


tfl  = 


L1  tan  o.{  tan2fl  +  tan4ftV/2(  tan2fl  -  tan2aV  /2(Z)  *  A") 

tan26> 


(4.27) 


al. 


cos 


'    \3/2 


-- r£— (i f^) 


sin  0 


a  j/1 


s  a 


(4.28) 


2,    ■  -1/       tan  a       »  ,        tan  a 
a2  =  a  {  sin    ( — r  )  + 


tan  8asp  |  |  tan  0 


flip 


1__tan^)|   — 
tan  tffl5p 


(4.29) 


a3  =  lab  \  D  x  h 
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where  L  is  the  height  of  the  cone  and  .V  is  the  normal  vector  of  the  surface  Av   The  de- 
tailed derivation  of  the  RCS  equations  is  given  in  Appendix  B. 

From  the  above  results,  the  RCS  value  of  the  projectile  in-flight  can  be  pre- 
dicted. It  is  shown  in  Figure  27  and  Figure  28  on  page  41.  From  above  figures,  it  can 
be  seen  that  the  maximum  RCS  is  about  0.09  m2  and  the  minimum  RCS  is  about 
0.02  m2.  Since  the  aspect  angle  of  the  projectile  changes  very  fast  right  after  it  is  fired, 
the  RCS  value  varies  rapidly  at  the  beginning  of  the  trajectory. 


Figure  27.      The  In-Fligl  t  Projectile  Radar  Cross  Section  vs.  Aspect  Angle. 
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Figure  28.      In-Flight  Projectile  Radar  Cross  Section  vs.  Slant  Range. 
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D.     PROPAGATION  EFFECTS 

Propagation  of  radio  waves  through  the  Earth's  atmosphere  is  significantly  affected 
by  the  propagation  medium.  The  presence  of  the  Earth  complicates  the  situation  in 
three  ways:  (a)  the  surface  of  Earth  scatters  radiation  and  produces  an  interference  pat- 
tern, (b)  the  Earth  casts  a  shadow  and  gives  rise  to  diffraction  phenomena,  and  (c)  the 
Earth's  atmosphere  is  inhomogeneous  and  produces  refraction  effects  [Ref.  13].  It  is 
customary  to  group  all  these  factors  into  a  single  quantity  called  the  propagation  factor 
(F).  In  this  thesis,  the  Engineer's  Refractive  Effects  Prediction  system  (EREPS)  [Ref. 
14]  is  utilized  to  compute  the  propagation  factor  based  on  local  environment. 

EREPS  uses  ray  tracing  for  radio  wave  propagation  within  the  horizon.  It  includes 
the  bending  effect  of  the  refractive  index  variation  of  air  and  the  divergence  factor  when 
the  ray  is  reflected  from  the  Earth.  The  field  at  the  target  is  the  sum  of  the  direct  ray 
and  the  reflected  ray.  resulting  in  interference  peaks  and  nulls  which  become  significant 
when  the  tarcet  falls  near  the  horizon  ( Fieure  29). 
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Figure  29.      Reflection  and  Interference  from  a  Spherical  Earth. 
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For  the  over-the-horizon  region,  the  waveguide  mode  approach  is  adopted  by 
EREPS.  The  Earth's  atmosphere  is  modeled  as  a  leaky  waveguide.  The  Earth's  surface 
is  the  lower  boundary  of  the  guide.  Modes  of  high  attenuation  along  the  horizontal  di- 
rection are  called  leaky  modes.  Modes  of  low  attenuation  are  called  trapped  modes.  In 
x-band.  the  effect  of  the  evaporation  duct  is  significant.  Trapped  modes  occur  when  the 
evaporation  duct  height  is  large  or  at  least  comparable  to  the  radar  wavelength  (see 
Figure  30  on  page  44)  [Ref.  15]. 

1.  Computation 

A  typical  EREPS  output  is  given  in  Figure  30  on  page  44.  The  local  annual 
average  meteorological  conditions  for  propagation  are  available  [Ref  13].  The  following 
parameters  are  used  : 

1.  Evaporation  duct  height  :    13.1m 

2.  k  (effective  Earth  radius  factor)  :    1.45 

3.  Wind  speed  :   6  m's 

4.  Absolute  humidity  :    lOg  m2 

5.  Surface  refractivity  :    370. 

Note  that  the  projectile  is  at  a  particular  height  when  it  arrives  at  a  range.  This 
computation  has  to  be  repeated  for  each  point  along  a  projectile  trajectory.  This  is  done 
for  several  trajectories  at  various  ranges  (Figure  32  on  page  45  and  Appendix  C  ).  Also 
note  that,  in  Figure  32  on  page  45.  the  ground  range  is  given  a  negative  sign  because  it 
is  measured  from  the  point  of  impact  toward  the  radar. 

2.  Summary  of  Results 

a.  Signal  Strength  Enhancement 

From  Figure  32  on  page  45.  it  can  be  seen  that,  with  ducting,  which  is 
present  most  of  the  time  at  the  test  site,  the  return  signal  will  be  substantially  enhanced 
from  the  corresponding  free  space  value.  (Note  that  the  computation  is  carried  out  for 
a  target  at  1  m  height,  i.e.,  when  the  projectile  is  near  its  point  of  impact.  The  applicable 
range  for  the  projectile  of  interest  lies  between  27  to  28  km.)  Thus  ducting  is  helpful  in 
enhancing  the  quality  of  the  recorded  full  range  tracking  data. 

b.  Strong  Interference 

As  shown  in  Figure  32  on  page  45,  the  return  signal  will  exhibit  strong  in- 
terference when  the  projectile  reaches  about  250  m  within  the  point  of  impact.  The  more 
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than  ±  15  dB  fluctuation  makes  it  difficult  to  track  the  projectile  in  this  region  using  the 
simple  CW  system  proposed  in  this  thesis.  On  the  other  hand,  since  this  strongly  fluc- 
tuating region  is  within  a  few  hundred  meters  from  the  impact  point,  this  small  portion 
of  the  trajectory  can  be  estimated  with  good  accuracy.  Hence  it  is  desirable  to  select  a 
'stop-processing  point'  at  which  the  CW  real-time  radial  velocity  and  range  computa- 
tions will  be  halted.  The  radar  antenna  will  be  guided  by  the  predicted  trajectory,  and 
the  return  video  sienal  will  still  be  recorded. 
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Figure  30.      Waveguide  Mode  with  Ducting  Height  and  Antenna  Height. 
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Figure  32.      Propagation  Factors  vs  Ground  Range  from  Impact  Point. 
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3.     Stop-Processing  Point 

It  should  be  noted  that  EREPS  does  not  take  into  account  the  antenna  radi- 
ation pattern.  Since  the  antenna  to  be  used  have  a  beamwidth  of  1°  only,  interference 
between  the  direct  ray  and  the  reflected  ray  will  not  be  important  until  the  antenna  ele- 
vation angle  is  lowered  to  within  0.5°  above  the  horizontal  direction  (Figure  33  on  page 
46). 
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Figure  33.      Narrow  Beamnidth  Effect  on  the  Spherical  Earth  Surface. 

This  criterion  can  be  used  to  determine  the  stop-processing  point  as  follows:  when  the 
system  records  an  antenna  elevation  angle  falling  below  0.5°,  the  square-law  detector 
and  the  comparator  will  be  switched  off.  The  servo  will  steer  the  antenna  according  to 
the  predicted  trajectory  smoothly  toward  the  projectile  point  of  impact.    Characteristics 
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for  two  such  stop-processing  points  for  a  radar-to-gun  distance  of  500  m  and  two  dif- 
ferent radar  heights  of  10  m  and  30  m  respectively  are  displayed  in  Table  1  on  page  47. 
Note  that  these  point  are  close  to  where  the  propagation  factor  starts  to  fluctuate 
strongly  as  predicted  by  EREPS. 

Since  the  processing  of  real-time  data  is  stopped  before  contributions  from  the 
reflected  ray  become  important;  in  what  follows,  the  propagation  factor  will  be  consid- 
ered as  unity. 

Table   1.     STOP  PROCESSING  POINT. 


Antenna  height  fm) 

K) 

30 

Radar  to  gun  distance(m) 

500 

500 

Projectile  height  (m) 

308 

305 

Slant  range  (km) 

27.53 

27.53 

Impact  range  (km) 

27.69 

27.57 

Ground  range  to  impact  point  (m) 

158 

157 

Elevation  angle  (degree) 

0.49 

0.45 

Elevation  angle  at  impact  point  (degree) 

-0. 1 5 

0.19 

Elevation  angle  difference  (degree) 

-0.64 

-0.64 

Azimuth  angle  (degree) 

36.5 

36.5 

Azimuth  angle  at  impact  point*  degree) 

36.9 

36.9 

Azimuth  angle  dillerence  to  impact  point(degree) 

0.39 

0.39 

In-flight  time  (sec) 

86.48 

86.48 
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E.     POWER  REQUIREMENT 

From  equation  (4.1),  the  power  requirement  for  tracking  up  to  the  stop-processing 
point  is  shown  in  Figure  3-4.  In  order  to  reduce  the  required  transmitter  power  re- 
quirement, a  multiple  filter  design  is  proposed  for  this  system. 

The  multiple  filter  design  is  based  on  the  characteristics  that,  after  an  initial  rise,  the 
doppler  frequency  decreases  with  increasing  slant  range  (Figure  35  on  page  49).  By 
switching  to  video  filters  of  narrower  bandwidths  at  longer  ranges,  the  coherent  inte- 
gration time  of  the  radar  is  increased  hence  the  power  requirement  can  be  lowered. 
Combining  the  fixed  range  independent  parameters  in  the  radar  range  equation,  the 
power  requirement  can  be  expressed  as  : 


(4.33) 


Pi 

=  (3.7339 

x  10  22)  x  ( 

R 

By 

a 

= 

cx(4- 

imaxJdmax  . 

c 

=  4.4S07  x  1CT22 

where  Bk  is  1. 

2/, 

(4.34) 
(4.35) 


Figure  34.      Power  Requirement  for  Full  Range  Tracking  without  Multiple  Filters. 
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Since  fd.  a  and  R  are  functions  of  the  trajectory  (Figure  35  on  page  49,  and  Figure  36 
on  page  50).  if  the  tracking  range  is  divided  into  three  regions  1,  2,  and  3,  and  one  filter 
is  used  for  each  region,  the  transmitter  power  requirement  becomes  : 


R 

't  ~  Q  ~~o~  1 1  ,max  ifd)\,max 

^V     O     lljmax  \Jd'2,max 
—  Q     a     ) 3, max  if dh, max 


(4.36) 
(4.37) 
(4.38) 


where  the  constant  C  is  the  same  as  that  of  equation  (4.34). 


Figure  35.      Doppler  Frequency  vs.  Slant  Range. 
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Figure  36.      Radar  Cross  Section  vs.  Slant  Range. 


The  power  requirement  written  in  terms  of  that  of  the  single  filter  case  is: 


Pt" 


a  -£— )     a 
c 


(4.39) 


where  C.  the  power  reduction  factor,  is  the  factor  by  which  the  power  requirement  is 
reduced. 

In  order  to  find  the  optimal  switch-over  points  for  the  filters,  the  logarithm  ofj^ 
versus  Rxfo,  normalized  to  (/^)m„  and  (IVIo)^  respectively,  is  plotted  (Figure  37  on  page 
51). 

Assume  that  the  first  switch-over  point  is  located  at  (xuyi),  and  the  second  one  is 
located  at  {x2,y2)-  These  values  can  be  determined  from  Figure  37  on  page  51  by  ob- 
serving that,  at  the  optimal  switch  points  : 


.r,  =  v,  +  -r-, 


(4.40) 


=  }'2 


(4.41) 
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The  optimal  locations  so  obtained  are  listed  in  Table  2  on  page  52.  The  power  re- 
duction factor  is  C.  =  10° "18  =  4.1S6,  as  seen  from  Figure  37  on  page  51. 

The  power  requirement  as  a  function  of  slant  range  when  three  filters  are  used  is 
plotted  in  Figure  38  on  page  52.  The  antenna  height  is  10  m  and  the  radar-to-gun  dis- 
tance is  500  m.  The  required  power  is  37  Watts.  A  conventional  CW  tracking  system 
will  track  to  the  summit  point  only.  This  point  is  located  at  a  range  of  region  1  and  re- 
quires a  transmitter  power  of  19.5  Watts.  Almost  full-range  tracking  is  possible  by 
doubling  the  conventional  power  through  the  use  of  multiple  filters. 
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Figure  38.      Power  Requirement  for  Trajectory  Tracking  with  Multiple  Filters. 


Table  2.     MULTIPLE    FILTERS    SWITCH    POINTS 

Parameters 

region  1 

region  2 

region  3 

stop  track- 
ing point 

Slant  range  (km) 

16.06 

23.73 

27.53 

27.53 

Maximum X  (kHz) 

S3  ?5 

22.03 

12. "2 

— 

Video  bandwidth 
(kHz) 

63.90 

26.44 

15.27 

... 

Radar  cross  section 
(m2) 

0.0385 

O.043S 

0.0862 

0.0903 

Required  power 
(Watts) 

36.29 

36.23 

36.29 

... 
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V.     CONCLUSIONS 

It  is  demonstrated  in  this  thesis  that  full-range  projectile  tracking,  with  real-time 
track  processing  stopped  at  a  range  within  160  m  of  the  projectile  point  of  impact,  can 
be  achieved  if  the  transmitter  output  power  is  37  watts  or  higher  and  a  set  of  three 
switchable  filters  are  used.  If  a  single  filter  is  used,  an  output  power  of  160  watts  is  re- 
quired. Either  configuration  can  be  implemented  with  a  commercially  available,  small 
size,  air-cooled    RF  amplifier  which  can  be  mounted  near  the  transmitting  antenna. 

The  tracking  relies  on  the  predicated  trajectory.  Although  real-time  dopplcr  infor- 
mation is  available  from  the  radar,  it  is  used  only  to  confirm  that  the  projectile  is  on 
track.  A  algorithm  to  utilize  this  real-time  information  for  updating  the  predicted  tra- 
jectory is  desired.  Kalman  filtering  appears  to  be  a  promising  technique  which  deserves 
further  investigation. 
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APPENDIX  A.  PROJECTILE  TRAJECTORY  SIMULATION  AND 

RANGE  ANALYSIS 


clear 
t2=clock; 
format  long  e; 

%  The  first  capital  letter  designations  are: 
%  G  for  the  gun  (muzzle);  R  for  the  radar. 


P  for  the  projectile; 


%   The  numerals  1,  2,    3  are  for  the  radial,  theta  and  phi  components  in 
%  an  inertial  spherical  coordinate  system;  4  is  for  the  radial  component 
%  in  the  cylindrical  coordinate  system  associated  with  the  spherical 
%  coordinate.   The  origin  is  at  the  center  of  the  Earth;  the  polar  (z) 
°0  axis  is  the  axis  of  spin  of  the  Earth. 

%  East  and  South  designations  are  local  quantities  in  a  frame  rotating 
%   with  the  Earth 


%  constant  declaration 
EarthT=86164. 06; 
EAV=2*pi/EarthT; 
a=( 63569 12+6378533) /2; 
g=9.81; 


%   sidereal  period  of  Earth  spin 
%   Earth  angular  velocity 
%   approximate  Earth  radius 
%  gravitational  acceleration 


%  latitude  24o36'38"  (radar  site) 
%   longitude  121o52'27"  (radar  site) 
?0  radar  transmitted  frequency  Hz 
%  radar  transmitted  wavelength  m 


RLATI=24+36/60+38/3600; 
RLONGI=121+52/60+27/3600; 

f=10525E6; 
lambda=3E8/f; 
?oinput  parameters 
! del  ww 
diary  ww 
RB=input( ' radar  to  gun  base  distance(m)  =?     ') 
%  RBaz=input(  '  radar  to  gun  direction  from  east=?    '),vpi/180 

RBaz=0; 
%   MN=input( 'Measurement  base  number  setting  =?     ') 
MN=1000; 

Rh=input( 'Radar  height(m)=?     ') 
%   Gel=input(  '  firing  elevation  angle(degree)=?     ')/180,vpi 

Gel=45/180*pi; 
%  Gaz=input( ' firing  azimuth  angle(degree)=?     '  )/180*pi 
Gaz=45/180--pi; 
diary 

answer=input( ' If  you  want  reenter  parameter  please  enter  1  else  enter'); 
if  answer=l 
return 
xx  a 
end 
Bl=5;  %   gun  barrel  length 


A=pi*(0. 155/2 )**2; 
m=45.  4; 
FAC=0.  5*A/m; 


%   projectile  cross  section  area 

%   projectile  mass  (kg) 

%   coefficient  for  computing  drag 
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%   measurement  base 

%   the  inertial  z-x  plane 


%  z-coordinate  of  radar 
%   radar  velocity 


%   muzzle  speed  of  projectile 


MB=MN*( lambda/2); 

%   Radar  Location  at  t=0 
Rl=a+Rh; 

R2=(90-RLATI)*pi/180; 
R3=0; 

R2s=sin(R2); 
R2c=cos(R2); 
R4=Rl*R2s; 
Rz=Rl-'-R2c; 
RV=R4*EAV; 

%   Firing  Direction  and  Speed 
PVo=850; 
Gels=sin(Gel); 
Gelc=cos(Gel); 
Gazs=sin(Gaz) 
Gazc=cos(Gaz); 

%   Gun  (muzzle)  Location  at  t=0 
Gh=Bl*Gels; 
Gl=a+Gh; 
Gl=Bl*Gelc; 

RGs=-RB-'-sin(RBaz) -Gl-Gazs; 
RGe=RB»cos(RBaz)+Gl-'-Gazc; 
G2=R2+RGs/Gl; 
G4=Gl*sin(G2); 
G3=R3+RGe/G4; 
Gz=Gl'vcos(G2); 

%  Displacement  Transformation 
G2c=cos(G2); 
G2s=sin(G2); 
DTG=( 0 , GAZC , GAZS;  0 , -GAZS , GAZC;  1 , 0 , 0-*-G2C , 0 , G2S;  0 , 1 , 0;  -G2S , 0 , G2C) ; 

%  Projectile  Initial  Location 
P1=G1; 

Pln=Gl; 
P2n=G2 
P3n=G3 
P4n=G4 


%   muzzle  height 


%  distance,  south  of  radar 
%  distance,  east  of  radar 


Projectile  Initial  Velocity 
PVln=PVo*Gels; 
PVg=PVo--"fGelc; 
PVe=PVg-Gazc; 
PAV2n=-  ( PVg-'-Gazs  )  /Gl; 
PAV3n=PVe/G4+EAV; 


%  velocity  along  ground 
%  velocity,  east 
%  angular  velocity,  theta 
%   angular  velocity,  phi 


%  Drag  Force  Considerations 


%  Wind  Profile:   assumed  constant  in  this  program. 
%     The  geostrophic  wind  should  be  measured  and  the  Ekman  spiral 
%     model  can  be  adopted  for  a  better  approximation  if  the  true 


55 


%  wind  profile  is  not  available. 
%   W ind= input ( 'Enter  the  wind  speed  m/sec  at  muzzle  height') 

Wind=6;  %   wind  velovity  at  muzzle  height 

%   WindBearing=input( 'Enter  the  true  wind  bearing  in  degrees') 

WindBearing=135;  %   clockwise  from  north 

Waz=((90-WindBearing)+180)*pi/180; 

We=Wind-'rcos(Waz);  %   wind  velocity,  east 

Ws=-Wind-sin(Waz);  %  wind  velocity,  south 

%   Temperature  Profile  (assumed  linear) 
%  Th=input( 'Enter  the  height  (m)  at  which  temperature  is  measured') 

Th=10; 
%   TREF=INPUT(' ENTER  THE  MESURED  TEMPERATURE  IN  CELSIUS') 

Tref=15; 

TrefK=273. 16+Tref; 
Thl=Th+a; 

S  =-0.0065;  %   TEMPERATURE  LAPSE  RATE  K/M 

%   TP=TREFK+S*(P1-TH1);  %   AIR  TEMPERATURE  AT  THE  PROJECTILE 

%  AIR  DENSITY 

%  DH=INPUT( 'ENTER  THE  HEIGHT  (M)  AT  WHICH  AIR  DENSITY  IS  MEASURED') 

DH=GH; 
X   DA0=INPUT( 'ENTER  THE  MESURED  AIR  DENSITY  IN  KG/CUBIC  METER') 

Dao=l. 225;  %  air  density  at  sea  level 

TDao=TrefK+s-(Dh-Th); 

ss=287. 05;  %  gas  constant  N-m/kg-K 

KC=1. 4;  %  Cp/Cv  for  air 

RAg=-l-g/(s"'ss);  %  calculation  factor 

%  DA=DA0*(TP/TDA0)**RAG;  %   AIR  DENSITY  AT  PROJECTILE  HEIGHT 

%   Speed  of  Sound 

KR=KC-ss; 
%   Vah=sqrt(KR*TP);  %   speed  of  sound  at  PI 

%   curve  fitting  for  Drag  coefficient  experimental  data  input) 
x=(  0,0.  05,0.  15,0.  2,0.  23,0.  25,0.  3,0.  35,0.4,0.45,0.5,0.55,0.  6,.  .  . 

0.  65,0.  7,0.  75,0.  825,0.  875,0.  98,1.  0,1.  05,1.  1,1.  2,1.  3,1.  4,1.  5,.  .  . 

1.  75  ,2  ,2.  5  ,2.  9); 

y=(0.'  16,0.  16,0.  16,0.  16,0.  16,0.  16  ,0.  16  ,0.  16  ,0.  16  ,0.  16  ,0.  16  ,0.  16, .  .  . 
0.  16,0.  16,0.  16,0.  16,0.  162,0.  17,0.  205,0.  225,0.  27  ,  0.  282  , 0.  29  , .  .  . 
0.  2875,0.  285,0.  275  ,  0.  25  ,  0.  23  ,  0.  195  ,  0.  17); 

cds=spline(x,y);  %   spline  curve  fitting  function 

%   get  spline  coefficient 

%   Trajectory  Computation  Begins 

t(l)=0;  %  initial  time 

n=0; 

^)'A-.VV.--;.--V-VV.-':'.-A'VrVrVr'^':V-.Vy.-*';V-:VVr   main  loop  begin  *"''"VV.-V.-T'-VfVr*Vr^^*T,rVrVr7WrVr*Vr*-***VrVr*Vr* 

while  Pln>a,  n=n+l; 

%  Update 
Pl=Pln; 
P2=P2n; 
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P3=P3n; 

P4=P4n; 

PVl=PVln; 

PAV2=PAV2n; 

PAV3=PAV3n; 

PV2=PAV2*P1; 

PV3=PAV3*P4; 

PVe=(PAV3-EAV)*P4; 

Height(n)=Pl-a; 

velocity  transformation  matrix 

PR3=P3-EAV*t(n);  %  location  relative  to  radar 

PG3=PR3-G3;  %   location  relative  to  gun 

P2s=sin(P2); 

P2c=cos(P2); 

PR3s=sin(PR3) 

PR3c=cos(PR3) 

PG3s=sin(PG3) 

PG3c=cos(PG3) 

VTRll=P2s*PR3c; 

VTR21=P2s*PR3s; 

VTR31=P2c; 

VTR=( VTR1 1 , P2c*PR3c , -PR3s;  VTR2 1 , P2c*PR3s , PR3c;  VTR31 , -P2s , 0) ; 


VTGll=P2s*PG3c; 
VTG21=P2s*PG3s; 

VTG31=P2c; 


X   Relative  Location 

PRXYZ=(P1*-VTR11;  VTR21;  VTR31)-(R4;  0;  RZ) 


SLR=NORM(PRXYZ); 


SLANT  RANGE  FROM  RADAR 


SRR(n)=SLR; 

°/  Vr  iV  ie  tfc  Vr  Vr  Vr  Vc  Vr  -,V  Vr  iV  ?V  if  ~V  Vr  ?V  Vr  ?V  Vr  Vr  Vc  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  >V 

/o 

DR(: ,n)=PRxyz/SLR;  %  d 

Dz=DR(l,n)*R2s+DR(3,n)*R2c;  %  t 

Dn=DR(3,n)*R2s-DR(l,n)*R2c;  %  t 

De=DR(2,n);  %  t 

traz(n)=atan2(Dn,De)*180/pi;  %  a 

trel(n)=asin(Dz)-180/pi;  %  e 

PGXYZ=DTG*(P1*-VTG11;  VTG21;  VTG31 
SL=norm( PGxyz); 
SRange(n)=SL; 
GRange(n)=a*acos((Pl**2+Gl**2-SL 

PGx=PGxyz(l); 

PGz=PGxyz(3); 

FRange(n)=a*atan(PGx/(Gl+PGz)); 

CRange(n)=PGxyz(2); 


sV  Vr  Vr  -.V  Vr  Vr  -V  Vr  - V  Vr  Vr  Vr  -V  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr  Vr 

irection  relative  to  radar 
racking  direction,  elevation 
racking  direction,  north  comp. 
racking  direction,  east  component 
zimuthal  angle 
levation  angle 

VrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVrVr 

)-(G4;  0;  GZ)); 
%   slant  range 

2)/(2*Pl*Gl)); 

%   ground  range 


%  range  in  firing  plane 

%  cross  range  from  firing  plane 


%Relative  Velocitv 
PV=(PV1;  PV2;  PV3); 
PVR=VTR*PV-(0;  RV;  0); 
VR(: ,n)=PVR; 
VRr(n)=PVl; 


%   velocity  relative  to  radar 


%   Time  Interval  and  Doppler  effect  Estimation 
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VDopp=PVR ' *PRxyz/SLR; 
Vdopl(n)=VDopp; 
fdl(n)=2*VDopp/ lambda; 
dtn=MB/VDopp; 
t(n+l)=t(n)+dtn; 
%   Frictional  Acceleration 


%   Doppler  velocity 

%   Doppler  Frequency 
%  square-law  doppler  detector 


TP=TrefK+s*(Pl-Thl); 

D  a=D  a  o  ••  ( TP  /  TD  a  o )  **RAg; 

Vah=sqrt(KR-''TP); 

PVW=( PV 1 ;  PV2 -WS; PVE -WE ) ; 

PVWM=norm(PVW); 

PWMM=PVWM/Vah; 

PH=PVW/PVWM; 

HeadingR(:  ,n)=VTR*PH; 

CD=ppval(cds,PWMM); 

FA=-FAC*Da*CD*PVWM*PVW; 

%   Location  and  Velocity  Prediction 


%   air  temperature  at  the  projectile 
%  air  density  at  projectile  height 
%   speed  of  sound  at  Pi 
%   PROJECTILE  AIR  VELOCITY 


%  projectile  heading 

%  heading  as  seen  from  radar 


%   drag  acceleration 


2; 


PA1=( PV2**2+PV3**2) /P1+FA( 1 ) -g*( a/Pl ) 

PAA2=P2s*P2c*PAV3**2+FA(2)/Pl; 

PAA3=FA(3)/P4; 

PVln=PAl*dtn+PVl; 

Pln=0.  5*(PVl+PVln)*dtn+Pl; 

PAV2n=(PAA2*dtn+PAV2)*(Pl/Pln)**2; 

P2n=0.5*(PAV2+PAV2n)*dtn+P2; 

P4n=Pln--sin(P2n); 

PAV3n=(PAA3*dtn+PAV3)*(P4/P4n)**2; 

P3n=0.  5*(PAV3+PAV3n)*dtn+P3; 
end 

%*************   end  of  main  loop  *************** 
%   Estimate  Point  of  Impact 

Ih=Pl-a; 

IA=-PA1; 

IV=-PV1; 

dt=(sqrt(IV**2+2*Ih*IA)-IV)/IA; 

N=n+1; 

t(N)=t(n)+dt; 

Pln=a; 

PVln=PAl*dt+PVl; 

PAV2n=(PAA2*dt+PAV2)*(Pl/Pln)**2; 

P2n=0. 5»(PAV2+PAV2n)*dt+P2; 

P4n=Pln*sin(P2n); 

PAV3n=(PAA3*dt+PAV3)*(P4/P4n)**2; 

P3N=0. 5*( PAV3+PAV3N)*DT+P3; 


%   Update 
Pl=Pln; 
P2=P2n; 
P3=P3n; 
P4=P4n; 
PVl=PVln; 
PAV2=PAV2n; 
PAV3=PAV3n; 
PV2=PAV2*P1; 


58 


PV3=PAV3*P4; 

Height(N)=0; 

%  velocity  transformation  matrix 

PR3=P3-EAV*t(N);  %   location  relative  to  radar 

PG3=PR3-G3;  %   location  relative  to  gun 

P2s=sin(P2); 

P2c=cos(P2); 

PR3s=sin(PR3); 

PR3c=cos(PR3); 

PG3s=sin(PG3); 

PG3c=cos(PG3); 

VTRll=P2s*PR3c; 

VTR21=P2s*PR3s; 

VTR31=P2c; 

VTR=( VTR11 ,P2C*PR3C , -PR3S;  VTR21 ,P2C*PR3S ,PR3C;  VTR31 , -P2S ,0); 

VTGll=P2s*PG3c; 
VTG21=P2s*PG3s; 

VTG31=P2c; 

%   Relative  Location 

PRXYZ=(P1*-VTR11;  VTR21;  VTR31)-(-R4;  0;  RZ)-; 

SRR( N)=norm( PRxyz) ; 
%********************  Radar  parameters  registration  ****************** 

DR(:  ,N)=PRxyz/SRR(N);  %   location  relative  to  radar 

Dz=DR(l,N)*R2s+DR(3,N)*R2c;     %   tracking  direction,  elevation 
DN=DR(3,N)*R2S-DR(1,N)*R2C;     %   TRACKING  DIRECTION,  NORTH  COMP. 
De=DR(2,N);  %  tracking  direction,  east  component 

traz(N)=atan2(Dn,De)*180/pi;    %   azimuthal  angle 
trel(N)=asin(Dz)"180/pi;        %  elevation  angle 

o/  y-  »•-  -*-  -♦.  Jm  «•-  *\.  -•-  -i- »'-  y.  .1- .'-  -»-  -t.  j-  -t-  -*  -  J-  .*-  .*.  j-  y- .  •-  y,  ■,%  -»-  -«-  j-  j-  -*-  j-  *»-  -•- »»-  -*-  *'r  *V  **■  *  V  V*  Vr  ~*r  ■*■  V*  ■*"  * V  *V  ***  *V  ^r  *V  *V  Vr  Vr  * V  Vc  Vf  Vc  ~V  *V  V?  tV  ~V  **■  V*  *V  V?  V  ?  V?  *V  Vr 

PGxyz=DTG*(Pl*[VTGll;    VTG21;    VTG31]-[G4;    0;    Gz]  ); 

SL=norm(PGxyz);  %  slant   range 

SRange(N)=SL; 

GRange(N)=a*acos((Pl**2+Gl**2-SL     2)/(2*Pl*Gl)); 

%   ground  range 
PGx=PGxvz(l); 
PGz=PGxyz(3); 

FRange(N)=a"vatan(PGx/(Gl+PGz) );     %  ground  range  in  firing  plane 
CRange(N)=PGxyz(2);  %  cross  range  from  firing  plane 

%  Relative  Velocity 
PV=[PV1;  PV2;  PV3]  ; 
PVR=VTR*PV-[0;  RV;  0]  ; 

VDopp=PYR' -''PRxyz/ SLR;  %   Doppler  velocity 

Vdopl(N)=VDopp; 

fdl(N)=2"VDopp/lambda;  %   Doppler  Frequency 

VR(: ,N)=PVR; 
VRr(N)=PVl; 

PW=(PV1;PV2-WS;PVE-VE); 
HeadingR( :  ,N)=VTR*(  PVW/norm(  PW)  ) ; 
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for  i=l:N; 

VRM(i)=norm(VR(:  ,i)); 

VRrf ( i)=VRr( i)/VRM( i) ; 

VRt(i)=(l-VRrf(i))*VRM(i); 

beamD(i)=l*pi/180*SRR(i); 
end 
EXETIMEl=etime( clock, t2) 
%****************  out  put  data  **^Wr**^Wf^Wr****************** 

diary  ww 
EXETIME1 

MN=MN 

ng=N 

MB=MB 

tg=t(N) 

r=Pl 

dlatitude=(P2-G2)*180/pi 

dlongitude=PG3*180/pi 

Rg=GRange(N) 

FRg=FRange(N) 

CRg=CRange(N) 

CRgmax=max( CRange ) 

maxHeight=max(  Height) 

maxfd=max( fdl) 

maxrange=SRR(N) 
VDOPmin=min( Vdopl ) 
VDOPmax=max( Vdopl ) 
diary 
**Vr*V:V:Vr******************  SAVE  AND  LOAD  DATA  FOR  NEXT  CALCULATION**  —  * 
!DEL  50010. MAT 

SAVE  50010  MAXRANGE  EXETIME1  DR  HEADINGR  SRR  HEIGHT  T  N  FD1  MB 
!DEL  X50010.MAT 
SAVE   X50010  TRAZ  TREL  VD0P1  VRM  VRRF  VRT  VR  VRR  T  SRR  BEAMD 

%-Vr*****-*-*****-********   Execute  Other  program  ********************** 

%PR  %   PRINT  OUT  TRAJECTORY 

%PRTRAK  %  TRACKING  PARAMETERS  &  OUTPUT 

%XX0  %  RCS  &  RANGE  ANALYSIS  PROGRAM 

%RCSA  %   RCS  SUBFUNCTION 

%RANGE3  %  POWER  REQUIREMENT 

OTATA  %  PROPAGATION  FACTOR  CURVE  FITTING 

0/  -U  »'-  - '-  -'-  -*-  -J.  J-  -U  »T-  .*-  JL  J*  -J.  JL  -•-  -'-  -•-  JU  -'-  -J-  J.  »'-  -'-  J-  J-  -»-  -*-  -J.  J-  ^  JU  -'-  »*-  .'-  -*-  -'-  -V  »'-  -t-  -'-  w3m  -'-  »*-  »'-  »'-  J-  .'-  -'-  ..'»  »'.  .'-  .*-  -'-  -'-  JL  -'-  -♦-  -»-  J-  -'-  -'-  -'-  -1-  .'-  -'-  ->-  -U  J-  ^ 


«£***************************outpUt   of   t r a j ectory"' 

! del   GH. met 

plot(GRange , Height) ,grid; 

tit le( 'PROJECTILE  TRAJECTORY'); 

xlabel( ' ground  range  (m)'); 

ylabel( 'projectile  height  (m)'); 

met a  GH 
%prtsc 
%pause 

! del  GC. met 

plot(GRange, CRange) ,grid; 

title( 'PROJECTILE  TRAJECTORY'); 

xlabel( 'ground  range  (m)'); 

ylabel( ' Cross  range  (m)'); 
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-  j*j-»'— »--,-j--,--j-.,-J-J-J-.*- 


meta  GC 
%   pause 
%   prtsc 

!  del  RH.  met 

plot(SRR, Height, V ) ,grid; 

title( 'PROJECTILE  TRAJECTORY'); 

xlabel( ' Slant  Range  (m)'); 

ylabel( 'projectile  height  (m)'); 
meta  RH 
%  prtsc 
pause 
%   !del  TH.  met 

plot(t, Height) , grid; 

xlabel( ' flight  time  (sec)'); 

ylabel( 'projectile  height  (m)'); 
meta  TH 
%  pause 
/oprtsc 

!del  TSGFC.met 

plot(t ,SRange,  t,GRange,  t,FRange,  t ,CRange) ,grid; 

xlabel( ' flight  time  (sec)'); 

ylabel( ' slant ,  ground,  firing  plane  &  cross  ranges  (m)'); 

meta  TSGFC 
%pause 
prtsc 

!del  SGFCT.  met 

plot( SRange , t ,  GRange,t,  FRange,t,  CRange ,t) ,grid; 

xlabel( ' slant ,  ground,  firing  plane  &  cross  ranges  (m)'); 

ylabel( ' flight  time  (sec)'); 

meta  SGFCT 
%pause 
%prtsc 

!del  TVR.  met 

plot(t,VRr,'+'),grid; 

xlabel( ' flight  time  (sec)'); 

ylabel( 'Vertical  Speed'); 

meta  TVR 
%pause 
%prtsc 
!del  TVH.met 

plot(t,VRt,'+'),grid; 

xlabel( 'Flight  time  (sec)'); 

ylabel( 'Horizontal  Speed'); 

meta  TVH 
%pause 
prtsc 
! del  FH. met 

p lot (FRange, Height) ,grid; 

xlabel( ' ground  range  in  firing  plane  (m)'); 

ylabel( 'projectile  height  (m)  ); 

meta  FH 
%pause 
%prtsc 
! del  FC. met 

plot(FRange, CRange) ,grid; 

xlabel( 'Firing  plane  Range  (m)'); 
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ylabel( ' Cross  range  (m)'); 
meta  FC 
%pause 
prtsc 

!  del  GV.  met 

plot(GRange,VRM) ,grid; 

xlabel( 'Ground  Range  (m)'); 

ylabel( ' Projectile  Speed  (m/sec)'); 

meta  GV 
%pause 
%prtsc 

! del  GRV. met 

plot(GRange,VRrf ) ,grid; 

xlabel( ' Ground  range  (m)'); 

ylabel( 'Ratio  of  radial  velocity  to  speed'); 

meta  GRV 
"opause 
%prtsc 

*****************  TRAKING  PARAMETERS  CALCULATION  AND  OUTPUT  ****** 
%****************  Tracking  angular  velocity  ************ 
avell=diff(trel). /diff(t); 
avazl=dif f(traz).  /diff(t); 

********************  PLOT  ****************************** 

for  n=l:  N-l 

ta(n)=t(n); 
end 

! del  AVAEt. met 

plot( ta.avell, '+' ) ,  grid; 

title( 'RADAR  TRACKING  ELEVATION  ANGULAR  SPEED  ' ); 

xlabel( ' In-flight  time  (sec)'); 

ylabel( 'Elevation  Angular  velocity  (degree/sec)'); 

meta  AVAEt 
%   pause 
prtsc 

! del  AVAZt. met 

plot( ta, avazl , '+' ) ,grid; 

title( 'RADAR  TRACKING  AZIMUTHAL  ANGULAR  SPEED'); 

xlabel( ' In-flight  time  (sec)'); 

ylabel( ' Azimuthal  angular  speed  (degree)'); 

meta  AVAZt 
°0pause 
prtsc 
%************************  Radar  Tracking  Parameter  output  ************ 

! del  tret. met 

plot(t,trel) ,grid; 

title( 'RADAR  TRACKING  ELEVATION  ANGLE  degree'); 

xlabel( ' In-flight  time  (sec)'); 

ylabel( ' Elevation  Angle  (degree)'); 

•meta  thel 
%   pause 
%  prtsc 

! del  trazl. met 

plot(t,traz) ,grid; 

title( 'RADAR  TRACKING  AZIMUTHAL  PARAMETER'); 

xlabel( ' In-flight  time  (sec)'); 

ylabel( 'Azimuthal  angle  (degree)'); 
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meta  trazl 
%pause 
%   prtsc 

!del  AVAE.met 

plot(avell, '+' ) ,grid; 

t it le( 'RADAR  TRACKING  ELEVATION  ANGULAR  VELOCITY  ' ); 

xlabel( 'MEASUREMENT  BASE  '); 

ylabel( 'Elevation  Angular  velocity  (degree/sec)'); 

meta  AVAE 
%  pause 
%   prtsc 

!  del  AVAZ.met 

plot(avazl, '+'  )  , grid; 

title( 'RADAR  TRACKING  AZIMUTHAL  ANGULAR  VELOCITY'); 

xlabel( 'MEASUREMENT  BASE'); 

ylabel( 'AZIMUTH  ANGULAR  (degree)'); 

meta  AVAZ 
%pause 
%prtsc 
EXETIMEl=etime( clock, t2) 

"/'jV^VVrVr's'riV-.VVrV.-Vr'jV-.'r-.'rVr'sVtV^V-.VVf-.'r-.'rVrV-    graph  'ViVVr^V-.ViV^V^^V^ViViVVrVr^ViVVriV'jV'^r'sV'jV'jViV'sV'sV'jV'sV'jV'jV'jV'jViVVr'jV'jV'jV 

!del  TVD.met 

plot(t  fdl),grid; 

title(fDOPPLER  FREQUENCY'); 

xlabel( 'Flight  time  (sec)'); 

ylabel( 'Doppler  Frequency  (Hz)') 

meta  TVD 
%   pause 
prtsc 

!del  treVD.  met 

plot(trel ,Vdopl) ,grid; 

title('DOPPLER  Velocity'); 

xlabel( ' Radar  tracking  elevation  angle  (degree)'); 

ylabel( 'Doppler  Velocity(m) ' ) 

meta  trehVD 
%pause 
%  prtsc 

! del  trazVD.  met 

plot( traz, Vdopl) ,grid; 

title( 'DOPPLER  Velocity'); 

xlabel( 'Radar  Tracking  azimuthal  Angle  (degree)'); 

ylabel( 'Doppler  Velocity(m) ' ) 

meta  trazVD 
%pause 
%prtsc 

!  del  RVDl.met 

plot(SRR, Vdopl) ,grid; 

title( 'DOPPLER  Velocity'); 

xlabel( ' Slant  Range(m)'); 

ylabel( 'Doppler  Velocity(m) ' ) 

meta  RVD1 
%pause 
prtsc 

!  del  tVDl.  met 

plot(t , Vdopl) ,grid; 

title( 'DOPPLER  Velocity'); 
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xlabel( ' In-flight  time(sec)') 

ylabel( 'Doppler  Velocity(m) ' ) 

met a  tVDl 
prtsc 
%************************  Radar  Tracking  Parameter  output  ************ 

! del  tret. met 

p lot (t.trel), grid; 

title( 'RADAR  TRACKING  ELEVATION  ANGLE  degree'); 

xlabel( ' In-flight  time  (sec)1); 

ylabel( ' Elevation  Angle  (degree)'); 

meta  thel 
%   pause 
prtsc 

! del  trazl. met 

plot(t ,traz) ,grid; 

title( 'RADAR  TRACKING  AZIMUTHAL  PARAMETER'); 

xlabel( ' In-flight  time  (sec)'); 

ylabel( ' Azimuthal  angle  (degree)'); 

meta  trazl 
%pause 
prtsc 
*  *  *  *  -.>  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

%*******  RCS  and  Power  requirement  Estimation******* 
load  50010 

t3=clock 

Evl=l. 2*max( fdl);  %  video  Bandwidth  for  1st  filter 

Bv2=l. 2*22031;  %  video  Bandwidth  for  2nd  filter 

Bv3=l. 2*12721;  %   video  Bandwidth  for  3rd  filter 

%****************  radar  parameters  calculation  ************************ 

Loss=10**0.  24;  °0  total  loss  2.  4dB 

KT=4 . 0  2 1 E - 2 1 ;  %  KT=K*T 

SN=10;  %  minimum  signal  to  noise  ratio  lOdB 

Lamda=3E8/10525E6;  %  wavelength  meter   10525MHz 

Fn=10**0. 2;  %   noise  figure   2dB 

G=io**4.43;  %  TX  or  RX  antenna  gain  44. 3dB 

const=( (4*pi)**3*KT*Fn*Loss*SN)/( (G*Lamda)2); 

O.   -'-  -'-  -•-  .'..  -'-  -'-  -*-  A  .*.  - '  -  mim  JL  J.  JU  ^.  *'-  -.'-  .'-  J-  J-  Jtm  -*-  *'-  JU  mtm  -'-  -'-  .'-  -'-  ~'-  .*.  -'-  **-  JL  -*-  -'-  -'-  -'-  -*;•  -'-  -'-  -'-  J-  ■.*-  -'-  -'-  -'-  J-  J-  J>-  Jm  -*-  J»  -J-  -A.  -'-  J -  ^-  J-  J-  J-  .'-  J-  J-  -*  -  JL  Jf  -I-  *ff  y- 

°b  Radar  cross  section  &  power  requirement  estimations 
for  n=l:  N 

if  max(Height)=Height(n)  %  Summit  point  data 

ttop=t(n);  %  time  at  summit 

Ltop=SRR(n);  %   range  at  summit 

end 
end 

^ *"'•  switch  1  power  requirement  '■•■* ««....■>. 

Bcon=Bvl*const; 
for  n-1:  1092 

SRRa(n)=SRR(n); 
ta(n)=t(n); 
uDl=DR(: ,n); 
h=HeadingR( : ,n); 
ct=h'*uDl; 

rcs=rcsa( ct);  %  res  function 

ctl(n)=ct; 
rcsl(n)=rcs; 
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%   power  requied  in  free  space 
%  required  power  at  summit 


%  res  function 


%   power  requied  in  free  space 
%  required  power  at  summit 


JU  JU  y-  A  -'-  J-  JU  -•  -  JU  -'r  JU  -'-  ~'  -  j-  -»-  J-  ^f  J-  J-  J-  J-  -'  -  J-  -'r 


thetal(n)=acos(ct)*180/pi; 

Ll=SRR(n); 

PtO=Ll**4*Bcon/rcs; 

Ptl(n)=PtO; 

if  Ll=Ltop 

ptop=PtO; 
end 
«£*******************  switch  2  power  requirement  **************** 

Bcon=Bv2"'vconst; 
for  n=1093:  1630 

SRRa(n)=SRR(n); 

ta(n)=t(n); 

uDl=DR(:  ,n); 

h=HeadingR( : ,n); 

ct=h'*uDl; 

rcs=rcsa( ct ); 

ctl(n)=ct; 
rcsl(n)=rcs; 

thetal(n)=acos(ct)*180/pi; 
Ll=SRR(n); 
PtO=Ll"-'-'4,-'Bcon/rcs; 
Ptl(n)=PtO; 
if  Ll=Ltop 

ptop=PtO; 
end 
switchl=max(  Pt 1 ) ; 
end  %   end  for  loop 

Bcon=Bv3*const; 
for  n=1631:  1898 

SRRa(n)=SRR(n); 

ta(n)=t(n); 

uDl=DR(:  ,n); 

h=HeadingR( :  ,n); 

ct=h'*uDl; 

rcs=rcsa(ct); 

ctl(n)=ct; 
rcsl(n)=rcs; 

thetal(n)=acos(ct)*180/pi; 
Ll=SRR(n); 
PtO=Ll**4*Bcon/rcs; 
Ptl(n)=PtO; 
if  Ll==Ltop 

ptop=PtO; 
end 
switch2=max(Ptl); 
end  %  end  for  loop 

«£**#*  ?rtr***********  ****  OU  t  put 

! del  wwl 

diary  wwl 

N=N 

f light time=max(  t) 

stoptime=max( ta) 

stoprange=max( SRRa) 

stopfd=fdl(1898) 

timetop=ttop 


%   res  function 


%  power  requied  in  free  space 
%   required  power  at  summit 


JL  JU  JU  JU  JU  J  -  JU  J-  -'-  JU  *'-  JU  -.' .  J*  -'-  J-  J -  J*  -'  -  J-  J-  J-  JU  -'-  J-  J-  J.  *  -  J-  -'-  J-  J:  J- 
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rangetop=Ltop 

powertop=ptop 

EXETIME2=etime( clock , t3 ) 

T0TALTIME=EXETIME1  +  EXETIME2 

diary 

%***********************  graph  radar  cross  section  ***************** 

! del  rcsll. met 

p lot (SRRa, rcsl, '+' ) ,grid; 

ylabel( ' Radar  cross  section  square  meters'); 

xlabel( ' Slant  range'); 

meta  rcsll 
%prtsc 
%   pause 
! del  rcsth. met 
plot(thetal , rcsl , '+' ) ,grid; 

ylabel( ' Radar  cross  section  square  meters'); 
xlabel( ' Aspect  angle  (degree)  ); 
meta  rcsth 
%   pause 
°0prtsc 

! del  rasp. met 
plot(ta, thetal , '+' ) ,grid; 

ylabel( ' Radar  cross  section  square  meters'); 
xlabel( 'Aspect  angle  (degree)  ); 
meta  rasp 
%pause 
?Dprtsc 

!  del  rasl.  met 
plot(ctl, '  +  ' ) ,grid; 
xlabel( ' in-flight  time  (sec)'); 
ylabel( ' Aspect  angle  (degree)'); 
meta  rasl 
°0prtsc 
%pause 

°  "-'""'"  **  *****  *****  ***  g  gy  g  variable  ******  ***************************** 

%!del  Pw50010.mat 

%save  Pw50010  ptlmax  ctl  t  fdl  rcsl  LI   N  thetal  Ptl  SRR  N 

O     Jm  mtm  .*..  mtm  mtm  m.*m  ml  m  Jm^m  m%m  m*m  »'j  mtm  »JL  mtm  »'■  mfm  mtm  mtmmUm  mfm  »>■  »  •  -l-L-l-t..l.,l,.».J„l,.l,J„l„>wl.1l,J..I„f, 

% r****c******    printout   power   requirement    •••••• * f 

! del  powerl.  met 

plot(SRRa,Ptl, '  +  '  ) ,grid; 

tit le( 'RADAR  POWER  REQUIREMENT  '); 

xlabel( ' Slant  range  (m)  '); 

ylabel( ' Power  requirement  (Watts)'); 

meta  powerl 
prtsc 
°0pause 

**********  RADAR  CROSS  SECTION  CALCULATION  ******************* 
clear 

LOAD  50010 
for  n=l:N 

uDl=DR(: ,n);  %   radar  to  target  unit  vector 

h=HeadingR(:  ,n);  /0unit  heading  vector  of  projectile 

ct=h',vuDl;  %   antenna  beam  angle 

%£ unction  [ res , ct] =RCS( ct) 
X   constant  declaration 

b=0. 155/2;  %  shell  base  radius 
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1=0.5236;  %  shell  cone(nuzzle)  height 

i=0.4144;  %  shell  cylinder  height 

ta=b/l;  %   =atan(b/l)  i.e.  half  shell  cone  angle 

ta2=ta*ta; 

ca2=l/(l+ta2); 

ca=sqrt(ca2); 

cx=abs(ct); 

cx2=cx*vcx; 

sx2=( l-cx2); 

sx=sqrt( sx2); 

rcsc=2,vi,vb"sx; 

b2=b*b; 

rcsb=pi,vb2,vcx; 
if  cx>=ca 

rcs=rcsc+rcsb; 
else 

rt=ta,vcx/sx; 

rt2=rt*rt; 

rcs2=b2-'-(asin(rt)+rt*sqrt(  l-rt2))*cx; 

rcsl=((b-->i)/(sx2)  )*(  l-cx2/ca2)**(  3/2); 

rcs=rcsl+rcs2+rcsc+rcsb/2; 
end  %  radar  cross  section  case  2 

rcs0(n)=rcs; 
ctl(n)=ct; 
thetal(n)=acos(ct)*180/pi; 

end 
! del  rcsll 

%   plot(SRR,rcsO, '+' ) ,grid; 

%ylabel( ' Radar  cross  section  square  meters'); 
%xlabel( ' Slant  range'); 
°imeta  rcsll 
°0prtsc 
%pause 
! del  rcsth 

plot( thetal , rcsO , '  +  ' ) ,grid; 

ylabel( ' Radar  cross  section  square  meters'); 
xlabel( ' Aspect  angle  (degree)  ); 
meta  rcsth 
%pause 
! del  rasp 

plot(ts thetal, '+' ) ,grid; 

ylabel( ' Radar  cross  section  square  meters'); 
xlabel( ' Aspect  angle  (degree)  ); 
meta  rasp 
%pause 
! del  rasl 

plot(ctl, '+' ),grid; 
xlabel( ' in-flight  time  (sec)'); 
ylabel( 'Aspect  angle  (degree)'); 
meta  rasl 
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APPENDIX  B.     CROSS-SECTION  ANALYSIS 

The  radar  cross-section  has  been  discussed  in  Chapter  4.  The  derivation  of  the  RCS 
equations  is  given  below.  The  radar  main  beam  direction  is  pointing  in  the  direction 
D.  The  projectile  axis  (or  heading)  is  pointing  along  h.  Accordingly,  the  coordinate 
system  (x,y,z)  can  be  chosen  so  that  x  =  D  and  the  y-z  plane  is  the  plane  on  which  the 
physical  cross-sectional  area  of  the  projectile  is  to  be  projected.  Another  coordinate 
system  (x',y'.z')  can  be  chosen  so  that  x'  =  h  and  the  projectile  is  a  body  of  revolution 
around  the  x  axis.   These  two  coordinate  systems  can  be  set  up  jointly  as  follows  : 

x  =  D  {B.\) 

:  =      Dx'h  {B2) 

\Dxh\ 

$  =  zxb  (B.y) 

x'  =  h  (B.A) 

z'  =  z  (B.5) 

y'  =  z'  xx'  (5.6) 

=  z  x  h  (B.l) 

[(/?'  x  D)h  -/)] 


Dxh 


(B.S) 


D  =      L       "  (B.9) 


rP  ~  r> 


Note  that  the  coordinate  (x.y.z)  are  used  in  this  Appendix  only  and  are  not  related 
to  other  parts  of  the  thesis. 

The  suface  of  the  projectile  is  given  by  : 
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y2  +  z'2  =fyx')  (5.10) 

The  RCS  is  the  area  enclosed  by  projecting  of  this  surface  on  the  y-z  plane. 

A.     CASE  1 

For  |  D  •  h  |  >  cos  a,  Figure  25  on  page  38,  the  conoid  is  hidded  in  the  shadow  re- 
gion.  The  projected  areas  can  be  easily  seen  to  be  : 

2 

Ax  =  3.14159-y-  (5.11) 

a,  =  A]  \D.h\  (BA2) 

A2  =  lab  (£.13) 

a2  =  A2\D*h\  (BAA) 

A2  =  A,  (BA5) 

o2  =  A3  \D*h\  (5.16) 


where  a  is  the  diameter  of  the  projectile  and  b  is  the  cylinder  height.  The  total  projection 
area  in  : 

o,  =  a,  +  a2  +  o3  (5.17) 

B.     CASE  2 

When  |  D  •  h  \  <  cos  a,  the  geometry  of  the  target  is  shown  in  Figure  26  on  page 
39.  The  RCS  value  consists  of  four  different  parts:  A}  ,  A2,  A3,  AA  .  Each  of  them  is  given 
below  : 

For  area  A{. 

y'2  =  cV2  (5.18) 

=  (x'-.x0')2ian7easp  (5.19) 

where 

c  =    tan  a  (5.20) 
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,/2    .    <\      t  ,„.  2n  i  ,2  « 2 


(c'  -  tan^flV).r"  +  2*0'  tan^flJ/rv'  -  V   tan^„.  =  0  (5.21) 


x0'  tan~0av, 
( tan2^5f  -  c2) 


.v,'c2  tan  Oasp 


where  xe',yc'  is  the  maximum  point  of  the  ellipse. 


y' 


(c  x'  +y'  tan  0ajp) 
(tan20a™  +  c4)1/2 


or 


Lc" 


tan  fl 


(5.22) 


JV  =  - — n it  (5-23) 

( tan  0a5_  -  c  ) 


(5.24) 


.v  tan  t>    r, 

J[x',y')  =  v'  tan  ^  +  x'  (5.25) 

A7  =  Vf  (5.26) 


(5.27) 


At  a  point  which  is  />,  ,  x'  =  L,  then: 

/  tan  Basp  +  Lc2  =  0  (5.2S) 


(5.291 


—  Lr2 
So  the  coordinate  at  n,  is  (-L  , : — ,  0  )  on  the  z   axis, 

tan  6csr 

if    2 

—Z-T +  z'2  =  c2L2  (5.30) 

tan  6asp 

or 

tan  flfl5;, 
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Hence 

Ax  =  zL{     X+C      )1/2  (5.32) 

tan% 

.  .,,  (tan20„D-c2)1/2 

=  L2c(  tan2^,  +  c4),/2 2f —  (5.33) 

tan  6asp 

Let  Ar  be  the  normal  unit  vector  of  surface  ^,: 

A  ->  A 

a  D  cos'a  —  /i  cos  6n,n 

.V=  f  (5.34) 

'   •    2Q  4,-4  2^, 

Y  sin  aasp  cos  a  +  sin  a  cos  #a5;, 

then  the  projection  area  for  ^,  is: 

L1  tan  a(  tan20a™  +  tan4a)1/2(  tan2fl^r  -  tan2a)1/2(/> .  N) 

ox  = P- - : (5.35) 

tan  6asp 

-> 
sm'6asp  cos  a 


A2  is  part  of  a  semi-circle  (Figure  26  on  page  39).   A  straight  forward  intergration  gives: 


A2  =  a2{s\n-\  ,    tan/    ,  )+  ,    ta"a    ;      /l ^^}  (5.37) 

2  l  I  ^n  %  |  |  tan  6asp  \  y]         tan^ 


2f    ■  -1/       tan  y       v  .        tana  /,         tan'a     ■,  .  A    l  i  /D1C1 


tan  #flJf  |  ;       |  tan  6asp 


tan  0dj;, 


^3  is  a  rectangle: 

/43  =  2a/>  (5.39) 


so 


a3  =  A3\Dxh\  (5.40) 


71 


where  a  is  the  radius  of  the  projectile  base,  and  b  is  the  height  of  the  cylindrical  part  of 
projectile. 

At  is  a  semicircle,  hence  : 

2 

AA  =  3.14159-—-  (£.41) 


a4  =  At\D.h  |  (£.42) 

The  total  cross-section  o,  is  the  sum  of  upper  4  projection  areas. 

ol  =  o,  4-  o2  +  a3  +  a4  (#-^3) 
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APPENDIX  C.     PROPAGATION  FACTORS 


Projectile 

Location 

Antenna  Height 

Ground  Rang 

to 
impact  poin 

e   Height 

10m 

L5m 

20m 

30m 

50m 

t 

(m) 

( 

m) 

(dB) 

(dB) 

(dB) 

(dB) 

(dB) 

-482.0000 

932. 

4000 

0 

0 

0. 

6000 

-0. 

3000 

-0. 3000 

-382. 0000 

741. 

2000 

0. 

6000 

0. 

6000 

0. 

6000 

-0. 

3000 

-0.  3000 

-333. 0000 

647. 

9000 

-0. 

7000 

-0. 

4000 

0. 

4000 

-0. 

3000 

-0. 3000 

-234. 0000 

455. 

2000 

-0. 

3000 

-1. 

1000 

1. 

2000 

1. 

2000 

-0.  2000 

-182.0000 

355. 

7000 

-0. 

9000 

1. 

0000 

1. 

7000 

0. 

7000 

1.  7000 

-152. 0000 

294. 

8000 

1. 

9000 

-3. 

7000 

1. 

2000 

-2. 

6000 

1.  2000 

-148. 0000 

289. 

2000 

2. 

3000 

0. 

5000 

-3. 

5000 

1. 

5000 

-3.4000 

-145. 0000 

283. 

7000 

1. 

5000 

1. 

0000 

0. 

3000 

0 

0.4000 

-142. 0000 

278. 

1000 

-0. 

7000 

1. 

3000 

-1. 

9000 

0. 

3000 

-0. 7000 

-140. 0000 

272. 

6000 

-3. 

8000 

-2. 

4000 

-1. 

0000 

-1. 

9000 

-2. 1000 

-137. 0000 

267. 

0000 

-3. 

0000 

-3. 

0000 

-3. 

0000 

-3. 

2000 

-0. 9000 

-132. 0000 

261. 

5000 

0 

1. 

4000 

2. 

3000 

2. 

6000 

2. 3000 

-131.  0000 

256. 

0000 

2. 

3000 

2. 

9000 

1. 

5000 

-2. 

4000 

-0.  3000 

-129. 0000 

250. 

4000 

3. 

2000 

0. 

7000 

-4. 

1000 

1. 

3000 

-2.  1000 

-125. 0000 

244. 

8000 

2. 

2000 

-5. 

1000 

0. 

2000 

0. 

2000 

2.  2000 

-123. 0000 

239. 

3000 

-0. 

8000 

-1. 

6000 

3. 

1000 

-1. 

4000 

-0.  7000 

-120. 0000 

233. 

7000 

-5. 

0000 

2. 

9000 

-4. 

0000 

-5. 

0000 

1.  6000 

-117. 0000 

228. 

2000 

-4. 

2000 

2. 

9000 

-4. 

0000 

-5. 

0000 

1.  6000 

-114. 0000 

222. 

6000 

-0. 

1000 

-1. 

0000 

2. 

9000 

3. 

2000 

0.  7000 

-111.  0000 

217. 

1000 

3. 

0000 

-6. 

5000 

2. 

7000 

-0. 

3000 

-0.  3000 

-108. 0000 

211. 

5000 

3. 

7000 

0. 

7000 

-5. 

3000 

0. 

8000 

-0.  3000 

-105. 0000 

206. 

0000 

2. 

7000 

3. 

7000 

-0. 

9000 

1. 

8000 

2. 7000 

-103. 0000 

200. 

4000 

-0. 

3000 

2. 

7000 

3. 

8000 

-2. 

4000 

-2. 4000 

-100. 0000 

194. 

9000 

-6. 

0000 

-2. 

2000 

0. 

6000 

3. 

0000 

-0.5000 

-97. 0000 

189. 

3000 

-5. 

3000 

-5. 

3000 

-5. 

3000 

-5. 

0000 

3.  0000 

-94. 0000 

183. 

8000 

-14. 

7000 

-4. 

7000 

-12. 

5000 

3. 

0000 

-11.  0000 

-91.  0000 

178. 

2000 

3. 

4000 

4. 

4000 

3. 

2000 

1. 

4000 

-1.  7000 

-88. 0000 

172. 

6000 

2. 

2000 

4. 

2000 

-3. 

5000 

-0. 

9000 

3.  1000 

-85. 0000 

167. 

1000 

3. 

5000 

-5. 

5000 

-1. 

5000 

3. 

5000 

-0.5000 

-83. 0000 

161. 

5000 

0. 

5000 

-3. 

5000 

4. 

2000 

-7. 

5000 

-4. 8000 

-80. 0000 

156. 

0000 

-6. 

4000 

3. 

4000 

2. 

3000 

4. 

1000 

2. 8000 

-77.  0000 

150. 

4000 

4. 

3000 

7. 

1000 

-10. 

0000 

-3. 

0000 

3.  0000 

-74. 0000 

144. 

8000 

0. 

7000 

0. 

9000 

2. 

3000 

1. 

9000 

-6.  9000 

-71.  0000 

139. 

3000 

3. 

7000 

-9. 

9000 

4. 

4000 

3. 

0000 

3.  0000 

-68. 0000 

133. 

7000 

4. 

6000 

-0. 

4000 

-1. 

4000 

-4. 

4000 

3.  2000 

-66. 0000 

128. 

1000 

3. 

8000 

4. 

4000 

-3. 

3000 

4. 

4000 

-6.4000 

-63. 0000 

122. 

5000 

1. 

1000 

4. 

1000 

4. 

1000 

-4. 

7000 

1.  9000 

-60. 0000 

117. 

0000 

-6. 

7000 

0. 

1000 

3. 

6000 

2. 

5000 

4. 0000 

-57. 0000 

114. 

4000 

-13. 

4000 

-4. 

4000 

0. 

2000 

4. 

6000 

4. 6000 

-54. 0000 

105. 

9000 

-0. 

2000 

0. 

2000 

0. 

2000 

4. 

1000 

2. 7000 

-52.0000 

100. 

3000 

4. 

0000 

4. 

8000 

4. 

7000 

4. 

7000 

3.  1000 

-49. 0000 

94. 

7000 

5. 

0000 

4. 

0000 

2. 

4000 

-3. 

9000 

-6.  9000 

-46. 0000 

89. 

2000 

4. 

5000 

0. 

2000 

-14. 

4000 

2. 

5000 

3. 6000 

-43. 0000 

83. 

6000 

2. 

8000 

-14. 

5000 

1. 

5000 

3. 

5000 

1.5000 

-40.0000 

78. 

0000 

1. 

4000 

0. 

2000 

4. 

7000 

-10. 

2000 

-2. 2000 
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37. 

0000 

72. 

5000 

-12. 

4000 

4. 

2000 

16. 

0000 

4. 

6000 

4. 

6000 

34. 

0000 

66. 

9000 

-4. 

8000 

4. 

7000 

-12. 

6000 

0. 

9000 

-3. 

8000 

31. 

0000 

61. 

3000 

1. 

4000 

1. 

4000 

1. 

4000 

-3. 

1000 

2. 

9000 

29. 

0000 

55. 

7000 

3. 

7000 

-7. 

0000 

4. 

3000 

4. 

5000 

2. 

6000 

26. 

0000 

50. 

2000 

4. 

7000 

-4. 

6000 

2. 

3000 

-1. 

9000 

-3. 

3000 

23. 

0000 

44. 

6000 

3. 

5000 

2. 

3000 

-6. 

4000 

-1. 

4000 

4. 

5000 

20. 

0000 

39. 

0000 

2. 

7000 

3. 

7000 

-2. 

1000 

4. 

4000 

-11. 

1000 

17. 

0000 

33. 

4000 

1, 

8000 

2. 

4000 

2. 

4000 

0 

4. 

5000 

14. 

0000 

27. 

9000 

1. 

1000 

2. 

4000 

3. 

7000 

-1. 

7000 

-9. 

0000 

11. 

0000 

22. 

3000 

-0. 

1000 

0. 

7000 

2. 

2000 

3. 

6000 

-10. 

7000 

-9. 

0000 

16. 

7000 

1. 

0000 

-0. 

3000 

0. 

3000 

3. 

3000 

4. 

1000 

-6. 

0000 

11. 

1000 

-0. 

6000 

-0. 

9000 

-0. 

2000 

-0. 

2000 

4. 

1000 

-3. 

0000 

5. 

5000 

4. 

7000 

3. 

6000 

2. 

7000 

2. 

4000 

2. 

8000 

0 

0 

-2. 

6000 

-4. 

8000 

-5. 

0000 

-4. 

6000 

-3. 

2000 
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APPENDIX  D.     MISSED  DETECTION  AND  FALSE  ALARM 

clear 

MN=1000; 

Npeak=MN/2; 

dtmax=0.  08650; 

Bvl=63900; 

FNnoise=Bvl"dtraax/Npeak; 

SNR0=10.  0; 

step=0.  5; 

Nf=5; 

Mf=29; 

thshstep=0. 005; 

thsh0=0.  355; 

for  m=l:Mf 

thsh=m"thshstep+thshO; 

thresh(m)=thsh; 

FNEnoise(m)=FNnoise-'2*asin(sqrt(thsh)  )/pi; 

end 

x=SNR0-step; 

for  n=l:Nf 

x=x+step; 

SNR(n)=x; 

sgnr=10--'"»-(0.  l*x); 

%  transfer  the  signal-to-noise  parameter  to  be  used  in  fdetect.  m 

snr=sgnr 

! del  hm  tempsnr. mat 

save  hm  tempsnr  snr 

for  m=l:Mf 
thsh=thresh(m); 
pfan=exp( -thsh-sgnr); 
Pf a(n,m)=pf an; 

FNerror(n,m)=pf an-'FNEnoise(m); 
end 

Pmiss(n, l)=Pmissdct( sgnr , 0 , thresh(  1) ); 
for  m=2:Mf 

Pmissd=Pmissdct(sgnr ,thresh(m-l) , thresh(m) ); 
Pmiss(n,m)=Pmiss(n,m-l)+Pmissd; 
end 
end 

!  del  hm  pfaulse.mat 

save  hm  pfaulse  Nf  Mf  thresh  FNEnoise  SNR  Pmiss  FNerror  Pfa 

!  del  hm  pfaulse 

diary  hm  pfaulse 

disp( '  Threshold   Missed      Noise     Faulse' ) 

disp( '  Detection    Error      Alarm   ' ) 

for  n=l:Nf 

SN_ratio=SNR(n) 

(thresh;  Pmiss(n,:);  FNerror(n, :  );  Pfa(n,:))'; 
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end 

diary 

pause 

!del  hm  pfaulse.  met 

plot  (thresh,  Pmiss( 1 ,:),'+' ,  thresh,  FNerror( 1 , : ) , ' o' ) ,grid; 

title( ' Signal-to-Noise  Ratio  =  10.0  db' ) 

xlabel( 'Threshold-to-Signal  Power  Ratio'); 

ylabel( 'Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  Pmiss(2 ,:),'  +  ' ,  thresh,  FNerror(2 , :  ) , ' o' ) ,grid; 

title( ' Signal-to-Noise  Ratio  =  10.5  db') 

xlabel( 'Threshold-to-Signal  Power  Ratio'); 

ylabel( 'Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  Pmiss(  3 ,:),'  +  ' ,  thresh,  FNerror( 3 ,:  ) , ' o' ) ,grid; 

title( ' Signal-to-Noise  Ratio  =  11.0  db' ) 

xlabel( 'Threshold-to-Signal  Power  Ratio'); 

ylabel( 'Missed  Detection  (  +  )  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  Pmiss(4 ,:),'+' ,  thresh,  FNerror(4, : ) , ' o' ) ,grid; 

title( ' Signal-to-Noise  Ratio  =  11.5  db') 

xlabel( 'Threshold-to-Signal  Power  Ratio'); 

ylabel( ' Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

plot  (thresh,  PmissC 5 ,:),'+' ,  thresh,  FNerror(5 , : ) , ' o' ) ,grid; 

title( 'Signal-to-Noise  Ratio  =  12.0  db' ) 

xlabel( 'Threshold-to-Signal  Power  Ratio'); 

ylabel( 'Missed  Detection  (+)  &  Noise  Error  (o)') 

meta  hm  pfaulse 

pause 

function  pmsd  =  Pmissdct( snr ,  thsha,  thshb) 

%   The  parameter  snr  for  the  singal-to-noise  ratio  (NOT  in  db)  should 

%     have  been  saved  to  hm  tempsnr. mat. 

?0!del  hm  tempsnr  mat 

/osave  hm  tempsnr  snr 

pmsd=quad( ' fdetect ' , thsha, thshb, 0.  0001); 

-■-  -•-  JU  -'  -  -•-  r'.  JL.  JU  J-  JL  .'.  JL  J.  -'-  J-  J*  JL  -■  -  .'  -  -•  -  .'  -  .'-  -•-  .' .  JL  -•.  JL  JL  J.  JL  JL  JL  JL  JL  JL  JL  .  ■-  JL  JLJUJLJL  JL  J  -  y  -  -'-  J.  JL  JL  JL  J-  J-  -'-  J  -  -'-  -'-  -J-  J/.  Jf  JU  JL  JL  J/f  JL  JU  .J.  JU  JL  JL 

%  The  parameter  sgnr  has  to  be  specified  in  the  program. 

function  pdens  =  fdetect(v) 

load  hm  tempsnr 

z  iv=2  -  snr "'•' s  q  rt  ( v )  *  i ; 

LIoz=log(BESSELN(0,ziv)); 

Lpdens=LIoz+log( snr) -snr-(v+l); 
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pdens=exp( Lpdens ) ; 

%*****************  end  of  program  ********************************** 

%*****************  output  data  ************************************* 
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%*  tV  iWoWr  •*  ?Y  ^  *  iV  •&  Vr  iV  *  •&  Vr  ^V  * 

end  of  data 
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